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Abstract 



We elaborate on the recently proposed static brane world scenario, where 
the effective 4-D cosmological constant is exponentially small when parallel 
3-branes are far apart. We extend this result to a compactified model with 
two positive tension branes. Besides an exponentially small effective 4-D cos- 
mological constant, this model incorporates a Randall- Sundrum-like solution 
to the hierarchy problem. Furthermore, the exponential factors for the hi- 
erarchy problem and the cosmological constant problem obey an inequality 
that is satisfied in nature. This inequality implies that the cosmological con- 
stant problem can be explained if the hierarchy problem is understood. The 
basic idea generalizes to the multibrane world scenario. We discuss models 
with piecewise adjustable bulk cosmological constants (to be determined by 
the 5-dimensional Einstein equation), a key element of the scenario. We also 
discuss the global structure of this scenario and clarify the physical properties 
of the particle (Rindler) horizons that are present. Finally, we derive a 4-D 
effective theory in which all observers on all branes not separated by particle 
horizons measure the same Newton's constant and 4-D cosmological constant. 
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I. INTRODUCTION 



Recent observational data [|lj indicate that our universe has a positive cosmological con- 
stant which, compared to the Planck or the electroweak scale, is many orders of magnitude 
smaller than expected within the context of ordinary gravity and quantum field theory. 
This is the well-known cosmological constant problem Recently, two of us proposed 

a static brane world solution to this problem in which the cosmological constant becomes 
exponentially small compared to all other scales in the model [|J|. Here, by "static", we 
mean the situation where the branes are stationary relative to each other []. We still allow 
for expansion of the Universe in the other three spatial directions, but confine ourselves to 
situations in which the Hubble expansion rate H is time-independent. In more complete 
Big Bang cosmological models, this would correspond to late times, when the cosmological 
constant becomes dominant. 

It is generally believed that fine-tuning is necessary for a very small cosmological constant 
in 4-dimensional theories JJ. This leads one to search for a naturally small cosmological 
constant in higher dimensional theories. However, for a usual compactification of a higher 
dimensional theory to an effective 4-dimensional theory, one ends up with a normal 4- 
dimensional theory, and the fine-tuning problem generically reappears. This is the case for 
usual Kaluza-Klein (KK) compactification, and for the generic compactifications with large 
extra dimensions |J. The Randall-Sundrum (RS) model |||7| provides a hope of avoiding 
this pathology. The RS model considers 3-branes inside 5-dimensional spacetime, with Anti- 
deSitter (AdS) spaces in the bulk. Normal finite gravitational interaction is reproduced even 
if the 5th spatial dimension is not compactified. For a single 3-brane (or a stack of 3-branes) 
with brane tension (or vacuum energy density) o"o, the effective 4-D cosmological constant 
as seen by observers on the brane is taken to be zero, that is MM 



A e// oc (T - ^/6A/k 2 = (1) 

where k 2 is the 5-dimensional gravitational coupling, and —A is the negative bulk cosmo- 
logical constant 0. We see that the brane tension, which is at least partly determined by 
the dynamics of the standard model, is not required to vanish. The cancellation between 
the brane tension and the bulk cosmological constant is imposed in any brane world sce- 
nario. However, the non-linear relation (i.e., a versus a/A), as well as the lack of the need 
to compactify the 5th dimension, follow from the non-factorizable metric in the RS model. 
Since the 5th dimension may stay uncompactified, any 4-D description will be inadequate. 
If one now insists on a 4-D description for observers on the brane, one ends up with a novel 
4-D description via the AdS/CFT correspondence 19,10], which offers hope to void the ar- 



gument for the need of fine-tuning to obtain a very small 4-D cosmological constant. This 



More precisely, the bulk spacetime possesses a timelike Killing vector field which is not hyper- 
surface orthogonal, so that the spacetime is stationary but not static, and the branes are invariant 
under the diffeomorphism generated by this vector field. 

2 Throughout this paper, the cosmological constants will have units of energy density in the ap- 
propriate dimension 
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has motivated a number of attempts to find a solution to the cosmological constant problem 



along this direction |[L1| . Here we shall elaborate on the approach of Ref ||. In the two 
brane Randall-Sundrum model |J, the cancellation in (JT[) is a fine-tuning. In addition, the 
two branes are required to have equal and opposite tensions. These two fine-tunings are 
removed in Ref [|J. 

The two brane model in Ref || does not solve the hierarchy problem, which is presumably 
solved in the RS model. However, the two brane model of Ref M is easy to generalize to a 
multibrane model in which both the cosmological constant problem and the hierarchy prob- 
lem can be explained simultaneously. The generalization to N 3-branes is straightforward 
because an observer on any brane can only see properties around his/her own neighborhood. 
This allows us to analyze the various possibilities generally. 

The simplest realization that solves both fine-tuning problems is a two brane compactified 
model. The model has an exponentially small A e jj, as in Ref ||, and implements the 
Randall-Sundrum-like solution to the hierarchy problem, but with a positive tension for 
the visible brane. Consider two parallel 3-branes with tensions er > a% > sitting in 
a compactified 5th dimension, with radius L2/271-. The o"o (Planck, hidden) brane sits at 
Lq = and the o\ (TeV, visible) brane sits at L\. Since L 2 is identified with Lq, this means 
the branes are separated by L\ on one side and by L 2 — L\ on the other side. Without loss 
of generality, let L 2 — L\ > L\. We find that, for large L±, the warp factor for the hierarchy 
problem is 

m Higgs _ A ( L i) _ -^(oo-oijzi/a /«\ 
m 2 Planck ~ A(0) ~ { > 

and the effective 4-D cosmological constant is 

_ 2 °o(°o + o-i) c _ K 2 [(ff[)+(ri)(L2 _ Ll)+((To _ Ji)Ll1/6 ^ ^ 
(cr - 

Since o"o > o\ > 0, and L 2 — L\ > Li, we see that 

ln ( m Planck \ > ln ( m Planck \ ^ 

\ A eff J \ m 2 Higgs J 

as is the case in nature, where \ia(mp lanck / A e ff) ~ 2.3 x 122 and ln(mp iancA ./m^ ggs ) ~ 
2.3 x 32. To get a feeling for the magnitudes of the various quantities, we may choose 
L 2 — Li ~ 2Li, and a ~ 2ai. If the brane tensions are around the Planck scale, the 
exponents that describe nature follow if L\ ~ ^0rripj anck . In this model, the branes are 
stationary (that is, they satisfy the brane equations of motion) and both L\ and L 2 are stable. 
But we still need a dynamical reason why L\ (or L 2 ) is large. However, if we understand 
the hierarchy problem via some means (e.g., comparing the renormalization group flows of 
couplings (marginal operators) versus the mass terms (relevant operators) in 4-D quantum 
field theory), this brane world model provides an explanation of the cosmological constant 
problem. As we shall see, the above inequality @ is robust in a generic multibrane world 
where the hierarchy problem is solved. 

The exponential behaviors are related to the warp factor in the Randall-Sundrum sce- 
nario, but with two key differences: 
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(1) The RS model demands A e ff to be exactly zero. This requires a fine-tuning. Here we 
find that this fine-tuning is not necessary for an exponentially small A e ff, as long as 
the branes are relatively far apart. 

(2) The RS model also requires the fine-tuning in (H). Not surprisingly, although k 2 (Tq— 6 A 
in the above model [Q is not zero, it turns out to be exponentially small. A priori, this 
is still a fine-tuning. To avoid this fine-tuning, the bulk A is treated as an "integration 
constant" that is determined by the 5-D Einstein equation. There are a number of 
well-known realizations where the cosmological constant is an integration constant 
[p~^ |20[1 . Since they are crucial to the above model, we shall review some of them, 



and then generalize them so that the bulk cosmological constant is piecewise constant, 
and each piece is an independent integration constant to be determined by the 5-D 
Einstein equation. We note that this idea can just as easily be incorporated into the 
RS model. 

The idea of treating the cosmological constant as an "integration constant" has a long 
history ||. Classically, it is not determined and so the cosmological constant problem is 
not solved. In the brane world, it is the bulk cosmological constant that is treated as an 
"integration constant" . In the single brane model, this bulk cosmological constant and so the 
effective 4-D cosmological constant are undetermined, so it takes a fine-tuning to suppress 
the 4-D cosmological constant. When the metric is non-factorizable, the introduction of 
another brane provides additional constraints which fix the value(s) of the bulk cosmological 
constant and hence the effective 4-D cosmological constant as seen by observers on the brane. 
In this sense, the multibrane world scenario is key to the solution. In Ref |4[], the bulk 
cosmological constant is treated as an integration constant in a generic model-independent 
way. Here, more details are worked out for the introduction of 5-form field strengths and 
its generalization. In specific models, we have additional equations, namely, an equation of 
motion for each brane. A priori, these equations are non-trivial for branes that are charged 
under the 4-form potential. We show that the static solution in Ref Q automatically satisfies 
these brane equations of motion. 

If the number of parallel 3-branes is not too large, sitting in the uncompactified (or 
compactified with large radius) 5th dimension, then we have a low density of 3-branes 
and their mean separation will be large in general. In the compactified version, observers 
on each brane will see an exponentially small effective 4-D cosmological constant. In the 
uncompactified case, all but one brane will have this property. In this sense, the smallness 
of A e ff is quite generic. 

In many cases, the metric between branes vanishes somewhere. These metric zeros are 
particle horizons which we show are coordinate singularities analogous to those in Rindler 
space. We eliminate those singularities, and deduce the global structure of the maximal 
extensions of the bulk spacetimes. We show that bulk regions separated by particle horizons 
often lie in separate connected components of the maximally extended solution, and therefore 
cannot physically influence one another 

In the multibrane model, an issue arises concerning the proper method for determining 
the 4-D Newton's constant Gn- The two standard methods to determine Gn are either to 
solve for the trapped gravity mode and Green's function, or to calculate the Hubble constant 
and determine the coefficient of the matter density term ||. Naively, one would expect the 
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normalization of the trapped gravity mode, and hence its wavefunction at a particular brane, 
to depend on the total number and placement of all of the branes, and so should the value 
of Gn found from the trapped mode. By contrast, the value deduced from the cosmo logical 
expansion rate should be determined locally, by the brane tension and the bulk cosmological 
constants on either side of the brane, and, one would think, should not depend on how many 
other branes there are, or where they are located. The resolution is that the normalization 
of the trapped gravity mode (or any other mode) should be computed only over the region 
between particle horizons. This also implies that the integration over the 5th dimension in 
the 5-D action to obtain the low-energy effective 4-D action |5| should be carried out over 
the same region. This means that we cannot compare the mass scales between branes that 
are separated by particle horizons, which are normalized independently. Any solution to the 
mass hierarchy problem can be addressed only between branes that are not separated by 
particle horizons. 

Now we can see the generic origin of the inequality (|J). To solve the hierarchy problem, 
both the Planck and the visible branes must be inside the same particle horizons. The cos- 
mological constant is exponentially small roughly as a function of the distance of the Planck 
brane from the particle horizon while the hierarchy scale (warp) factor is exponentially small 
roughly as a function of the distance of the Planck brane from the visible brane. Since the 
visible brane must be between the Planck brane and the particle horizon, the inequality (|J) 
follows. 

Using the effective 4-D action approach, we show that both Gn and the effective cosmo- 
logical constant A e ff are the same for all observers on all branes within a pair of particle 
horizons, irrespective of the type of brane (positive or negative tension, Planck or visible) on 
which they live. Generically, the visible brane will have different bulk cosmological constants 
on its two sides. The correction to the Newton's law is calculated in this case; it has the 
same dependence on particle separation as in the symmetric case j7|, but with a modified 
coefficient. 

In general, A e //(Lj) is expected to be more complicated than given above, because we 
expect additional non-gravitational brane-brane interactions at small separations. For large 
separations, it is reasonable to assume that the inter-brane dynamics is dominated by pure 
gravity as described here. In a more realistic situation, the matter density on the visible 
brane (and dark matter on the other brane) should be included, and the branes should be 
allowed to move. We shall not consider these effects in this paper. However, we do note that 
if the kinds of explanations offered here are indeed behind the smallness of the cosmological 
constant and Higgs mass relative to the Planck scale, branes must have moved very little 
since the epoch of cosmological nucleosynthesis. Otherwise, the Higgs mass scale might have 
changed enormously, which would be inconsistent with the success of Big Bang cosmology 
in explaining the light element abundances. 

The paper is organized as follows. In §j[]j, we review models in which the bulk cosmo- 
logical constant emerges as an "integration constant" and discuss the generalization where 
it becomes piecewise constant and piecewise adjustable. In §|I1 we present the multibrane 
solution of the model and review the 4-D effective action. § |TV| reviews and elaborates on 
the two brane model of Ref [ffl. §|V| considers a two brane compactified model where the 
Randall- Sundrum solution to the hierarchy problem is incorporated together with an ex- 
ponentially small effective 4-D cosmological constant. § |VT| gives a general analysis of the 
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multibrane world. §|VT| discusses the global structure of the metric. In § |VIII| , we discuss 
the implications of the particle horizons on the determination of and the mass hierarchy 
issue. §|X| considers the correction to Newton's gravity law in the case where the branes 
are sitting between two different AdS spaces. §|X| contains some overall discussions and §|X| 
gives a brief summary. Some of the details are relegated to the appendices. Appendix |A] 
solves the brane equation of motion for a stationary brane. Appendix [B] derives some de- 
tails of the two brane compactified model. Appendix y gives the details of the coordinate 
transformation used in the analysis of the global structure. For the visible brane in Ref 
or the two brane compactified model, a naive determination of Newton's constant using the 
Hubble constant approach sometimes yields a wrong result. In Appendix |B], we comment 
on the ambiguity /problem with this naive approach. Appendix [E] gives a Feynman diagram 
analysis of Newton's force law on the relation between different metrics. 



II. BULK COSMOLOGICAL CONSTANTS AS INTEGRATION CONSTANTS 

One of the key points underlying the two brane model [Q is that the bulk cosmological 
constant is not an input parameter of the model, but an "integration constant" to be deter- 
mined by the 5-D Einstein equation for the bulk together with the boundary conditions at 
the branes. 

Suppose the 5-D theory discussed above arises from the compactification of a higher 
dimensional theory. Then the bulk cosmological constant, the brane tensions, as well as the 
brane charges depend on the various compactification radii. If we do not fix the compacti- 
fication volume but instead make the compactification radii dynamical, this will provide a 
way to adjust, among other quantities, the bulk cosmological constant. This means that the 
bulk cosmological constant becomes dynamical. We will treat possible dynamical models 
for the brane world elsewhere. 

Here we shall review two approaches in which this "integration constant" appears: 



i) unimodular gravity ||12|-[L4||, which is suitable for the two brane model; 



(ii) 5-form field strength model ||17|- f20fl , which is suitable for the multibrane model, either 
compactified or uncompactified, but not suitable for the orbifold version. 

We then present a generalized model. Special cases of this generalized model reduce to the 
5-form field strength model and the fully covariantized variation of the unimodular gravity 
16fl . To adapt it for the multibrane orbifold model, we introduce brane couplings so that the 



bulk cosmological constant becomes piecewise constant. We note that 5-form field strength 
also appears in gauged supergravity and/or superstring realizations of the RS scenario plf . 



In 4-D gravity, when the cosmological constant is an integration constant, it is typically 
left undetermined. In the multibrane model, the bulk 5-D cosmological constant is piecewise 
constant, determined by the 5-D Einstein equation, including the jump conditions at the 
branes and the brane equations of motion. As we shall see, the number of such adjustable 
bulk cosmological constants should be equal to the number of bulks between branes, so, for 
N branes in either the uncompactified case or the S 1 /Z 2 orbifold case, we require N — 1 
adjustable bulk cosmological constants. In the compactified version, we need N adjustable 
bulk cosmological constants. 
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A. Unimodular Gravity 



It was first pointed out by Einstein [12] that the cosmological constant could arise as 
an integration constant if the original equations of general relativity were replaced by their 
trace-free forms. Einstein [|L2j had in mind traceless source terms due to electromagnetic 
fields alone, but the same traceless field equations are also obtained in "unimodular gravity" 
13|,15l, where the conformal part of the metric is constrained and is not allowed to vary. The 



fields in this theory consist of a metric g a t, a fixed, background, non-dynamical volume form 
and other matter fields. The action of the theory is the standard Einstein action that 



-a\a2-..a d i 



depends only on the metric and the matter fields. However, only metrics g a b whose volume 
forms coincide with the background volume form e are allowed. (An equivalent description of 
the restriction on the space of metrics is that only metrics whose determinants in a particular 
fixed coordinate system are —1 are allowed |13[)- Unimodular gravity is consistent with a 
massless spin-2 graviton propagating with ((d— l){d— 2) — 2)/2 polarizations in the linearized 
version of gravity in ci-dimensions ||13|| . Varying the action over the allowed class of metrics 
yields in d dimensions only the traceless part of Einstein's equation, 



R 



ab 



g ab R 



d 



T 



ab 



(5) 



where T = T°. 

The standard argument of invariance of the action under arbitrary linearized coordinate 
transformations still applies and shows that energy-momentum tensor satisfies the usual 
conservation law, 



ab 



0. 



and of course the Bianchi identities, 



1 



V b ( R ab - -g ab R ) = 0, 



still hold. Combining Eqs. © - (0) gives 

V a \(d/2-l)R + k 2 T 



0. 



(6) 



(7) 



(8) 



Equation ([8]) implies that (d/2 — 1)R + n 2 T is a constant, which we will denote by n 2 Ad. 
As was noted by Einstein [0, when the relation (d/2 — 1)R + k 2 T = n 2 Ad is inserted in 
Eq. (H), the result is the full Einstein equation with cosmological constant A: 



G ab + K z Ag ab = k 2 T 



ab 



Jlrpab 



(9) 



In this model, the cosmological constant A is merely an integration constant, and has nothing 
to do with any input parameter in the action, the microphysics, or the (quantum) vacuum 
fluctuations. (Note that we may absorb any cosmological constant term inside T ab into 
A.) In a set-up where A is not determined, it becomes a free parameter. In standard 
4-D gravity, this altered status of the cosmological constant does not help to explain its 
smallness. However, in the two brane world scenario in 5-D, A is no longer a free parameter. 
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Rather, it is the bulk cosmological constant determined by the 5-D Einstein equation and 
the boundary (jump) conditions at the branes Q 

The presence of a non-dynamical background field in unimodular gravity is of course an 
unattractive feature of the theory. However, it is possible to find other theories with no 



background fields which, classically and on-shell, are equivalent to unimodular gravity flU 
For example, one can treat A as a dynamical field and introduce an appropriate Lagrangian 
multiplier (namely, a vector field) to render it spacetime- independent [IB|. This theory is a 
special case of our generalized model of Sec. [II C| below. 



B. 5-Form Field Strength 



Another way to get a piecewise constant bulk cosmological constant is to introduce 
(d — l)-form fields |T7^pO| in d-dimensional spacetime. This scenario is natural in compact- 
ified versions in string/M theory, in which D-branes are charged. This permits multibrane 
scenarios with the cosmological constant taking different constant values in different regions 
of the spacetime. The discontinuities are due to charged (d — 2) -branes, which act as sources 
for the (d — l)-form field. A worldvolume action for these branes should, therefore, contain 
a Wess-Zumino term: the integral of the pullback of the (d — l)-form potential. In the 
multibrane world, the action S consists of a bulk action containing gravitywith metric g a b, 
coupled to a (d — l)-form potential A^-i), with field strength = dA^-i), and a set of 
standard worldvolume actions of (d — 2)-branes containing the above-mentioned WZ terms 
with dynamical coordinate fields X® This means the location of the nth brane in the bulk 
is given by the embedding functions x a = X%(££). In addition, the nth brane possesses a 
brane metric j£ v which is a function of the brane coordinates ££.(for d — 5, a — 0, 1, 2, 3, 5 
and [A = 0, 1, 2, 3). The action for such a system is [p2^ , p0f 



S = J d d x^f\g] 



+ E { -ifn J d^^\ ln \ l^d,X a d u X b g ab (X n ) - C 



(d-1) 



j / S d - l ^ n A ai ... ad ^ (x n ) d^x* ■ ■ ■ ^.r-^ 1 -^ 1 1 



(10) 



where a n and e n are the tension and the charge of the nth brane. The dimensionless constant 
C must have the value p — 1 in order to recover the correct equations of motion (see below). 
The tensor density e^ 1 -^- 1 is totally antisymmetric with e 01 -- ( d - 2 ) = i. Varying the action 
S with respect to the metric, we have 



G ab = k 2 { -kg ab + 



ab 



1 



(d-1)! 



Cl - Cd - 2 ~ 2d 9 {d) 



E 



Or. 



d (d-1 Wl7n| l^d^X a d v X b 8^(x-X n 



(11) 



Varying the action 5" with respect to A^-i), 
we have, respectively, 



the brane metric and the coordinates X a , 



S 



;i2) 



9fj,X^d u X n g a b (X r , 
and the brane equation of motion for the nth brane 



(13) 



+ 



(d-r 



T?a fx ybi fx yb 



0. 



(14) 



Eq. (|T3|) simply states that the worldvolume metrics are those induced on the worldvolumes 



by the embedding coordinates X". In equation Eq. (13) there is a coefficient C / (p — 1), for 
a p-brane (p = 3 in our case), which we have set equal to 1, i.e., is the induced metric 
on the brane. Had we written the brane tension term in the Nambu-Goto action form, the 
same equations would have resulted. In the bulk 



d a [J\g\F 







which has solution 



(15) 



(16) 



where e is a constant that can vary from one interbrane region to another, and the tensor 
density e ai "' ad is totally antisymmetric with e 01 — = l. Inserting this solution in Eq. (|TTD 
we see that the only contribution of the field strength is to the bulk cosmological constant A. 
In this paper we are intersted in the parallel stationary 3-branes case, so we may choose the 
Gaussian normal coordinate: ds 2 = g^ v (x,y) dx^dx v + dy 2 . The brane equation of motion 
(n|) for a stationary brane is discussed in Appendix |A]. Using worldvolume reparametrization 
invariance we can set d—1 coordinates (static gauge) to the values X® = 5°££. The remaining 
coordinate, namely X d = y, is along the dth direction, 



xi = L r 



(17) 



where the L n 's are constants. We can perform the volume integrals in Eqs. ( |T^ , [T2| ) leaving 
only 1-dimensional delta functions 5{y — L n ). Also, this implies for the worldvolume metrics 
that: jn^u = g^ u (L n ). This can also be seen by solving the equation Eq. (|IB|) in the 
Gaussian normal system of coordinates in the bulk. Since the branes are parallel, the 
geodesies orthogonal to one brane will be orthogonal to all the branes. We obtain: 



7n/iK — d^dv^dab (X n ) — g^ (L n ) 



(18) 



where we have taken g§\ = for Gaussian normal coordinates. It follows that y |7| = y\g\, 
and ( |T2| ) becomes 



\ g \V a F°M~b^i _ £ e n e b ^ b ^5 (y - L n ) = 0. 



(19) 
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The determinant of the metric is a continuous function everywhere in the space, so as 
we go across the nth brane with charge e n , the field strength jumps by e n , from e( n ) to 
C( n +i) — e (n) + e n . In the multibrane world, the resulting cosmological constant in the bulk 
between the (n — l)th and the nth branes is given by 



- A n = A + -e( n) = A + - (e (0) + e + e x + ... + e n -\ 



(20) 

where e(o) is a constant background field strength. Now the Einstein equation ( |TTD becomes 

" N N-l 

Y: A n e(y - L n ^)e(L n - y ) g ab - ]T -fsffirs ( y - u 



G ab = k 2 



n=0 



n=0 



(21) 



where A n > for AdS space, with — L_i = = oo in the uncompactified case. The 
multibrane world with an uncompactified 5th dimension is shown in Figure 1. 
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FIG. 1 Schematic multibrane setup. 

Although the charges e n of the branes are taken to be fixed, the background field strength 
e(o) is an integration constant (to be determined by the Einstein equation.) To have AdS 
spaces in the bulk, we require A to be negative enough so that all A n > 0. Since the 
(d — l)-form field has no dynamical degrees of freedom, it is consistent to flip the F 2 term in 
the action S so that the field strength contributes only negatively to the bulk cosmological 
constant, avoiding the need to introduce a negative A. However, if the (d — l)-form field 
arises as a component of an antisymmetric field in higher dimensional space, then positivity 
of its higher dimensional kinetic energy density does not allow this freedom. 

Suppose the <ith dimension (i.e., the y direction) is compactified in the N 3-brane 
scenario. Boundary condition requires the sum of the brane charges to be zero, that is, 
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eo + ei + e<i + ... + = 0, independent of the value of e( ). If the model is a S 1 /z 2 orbifold, 
the symmetry constraint requires that = 0. In this case, the bulk A n are interrelated, 
so we have to find another way to allow A n to be adjustable. 

In the two brane world, e( ) is determined by the Einstein equation. In the multibrane 
world, we may need more freedom than one integration constant. To achieve this, we can 
consider a scenario where there is more than one (d — l)-form field. Considering the case 
d — 5, the introduction of more than one 4-form field is quite natural in string/M theory. To 
be specific, consider M theory which has 2-branes and 5-branes, which are electrically and 
magnetically charged, respectively, under the 3-form field A( 3 y The dual of this 3- form field 
is a 6-form field A^y If 6 of the 11 dimensions of M theory are compactified toroidally, then 
of the 6-form field there are, among other fields, 15 4-form fields Ah (J = 1, 2, 15) 
in the uncompactified 5-D spacetime. The 3-branes in this remaining 5-D spacetime are 
5-branes of the 11-D theory with two of their spatial dimensions compactified. There are 
15 types of such 3-branes, each charged under one 4-form field A J ^. In more realistic 
compactifications, some of these 4-form fields may be projected out, although generically, 
we expect a number to remain. 

Let each 3-brane be a stack composed of these compactified 5-branes. Suppose we have 

M 4-form fields Ah and their corresponding field strength Fh, J = 1,2, , M, under 

which the nth 3-brane has charge e^, J = 1,2, , M. Then the resulting cosmological 

constant in the bulk between the (n — l)th and the nth branes is given by 

- An = A + \ E«)) 2 = A + \ E «o) + e J +e{ + ... + e^) 2 . (22) 
j j 

Again we require A„ > for AdS spaces between branes. Now there are M background field 

strengths eh, J = l,2, , M, that are integration constants to be determined. Depending 

on the charges of the brane, this will allow some number of bulk cosmological constants to 
be adjusted to satisfy the Einstein equation. 

C. A Generalized Model 

We are interested in models where the bulk cosmological constant is piecewise constant, 
with each piece an integration constant that will be determined by the Einstein equation. 
The cosmological constant in the unimodular model is the same everywhere. This scenario is 
suitable for the two brane model, but is inadequate for the multibrane world when there are 
more than two branes. The model with M types of 5-form field strengths is suitable for the 
(M + l)-brane model, in the uncompactified case, or the M-brane model, in the compactified 
case. However, it is not suitable for the S 1 /z 2 orbifold model, since the background 5- 
form field strength must be set to zero due to the Z 2 symmetry. As a consequence, the 
bulk cosmological constants are fixed, so it is useful to consider scenarios where the bulk 
cosmological constant can be piecewise adjustable. Here we present such a model, which is 
a generalization of the above models. 

Let us first consider a model without branes, which includes a scalar field and its 
effective potential, A(0): 

S = / A V^{^- A ^)- r ^}' ( 23 ) 
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The independent variables of the model are g a b(%), T a (x) and 0(x); the Einstein equation 



is 



G afe = /€ 2 [-A(0)-T c 9 c 0]^ fe . (24) 
The variations of S with respect to T a and give the following equations: 

SS 



5T a 







g\d a 







(25) 



g-o^*(,/5P-)-,/w£ 

Eq. implies that (x) = constant, so that ( [Zip becomes 

G a b = —A (0 = constant) g a b 



(26) 



(27) 



where A(0) is an integration constant. Let us now consider special choices of A(0). For 
example, if 



A (0) = -0 2 , 



(28) 



the action S reduces to the <i-form field strength case we have already discussed |17| . To see 
this, define T a as the dual of a (d — 1) form A^-i)'- 



rpa 



c ab\---b d . 



(d-l)W\g 



A 



bi-b d 



(29) 



Then Eq. (|26| ) becomes 



c ab\---b d _± 



(d-iy. 



9 a -Afti— bd-j 



r ai— a d 



P 



(30) 



and substituting this back into the action (or equivalently, perform the functional integration 
over in the generating functional) and using (e ai ' ad F ai ...a d ) 2 = — (\g\)d\F ai ...a d F ai '" ad , gives 

ld,x M^ R -2^ F °-°< pn ""}- (31) 

Alternatively, if we choose 



A(0) = 



(32) 



we recover the model of Ref ||16|| , which may be considered as a covariantized version of the 
unimodular gravity. The value of can be arbitrary in the orbifold model. We may also 
combine the two approaches to obtain an action S 



(33) 



S = J d^x^l^R-^Fa^F^-^-T^a^ 
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We can easily further generalize this formalism to M fields Tj and <pj, J 
couple them to the (d — 2)-branes. 



1,2 ---M and 



S = 

E 



d {d) xy\g 



1 M 1 

—R-J2<pj + T]d a <Pj\ + 
jK j=i J 



Tn 



^ v d^X a d u X h g ab (X n ) - C 



i.J 



(34) 



If the nth 3-brane has charges u^, J = 1,2, , M, then the resulting cosmological constant 

in the bulk between the (n — l)th and the nth branes is given by 



A, 



E< 

j 



E 

j 



(35) 



Again we require A n > for AdS spaces between branes. There are M background constants 

<fr J , J = 1,2, , M, that are integration constants to be determined. In this case, only the 

sum of 0's is determined by the Einstein equation, so only one bulk cosmological constant 
can be adjusted to satisfy the Einstein equation. Here, A(0) = can be either positive or 
negative. We can combine one <fi field with M 4-form potentials. We shall comment on the 
dynamics of this system when applied to the static multibrane solution. 



III. SOLUTION TO EINSTEIN EQUATION AND GENERAL FORMALISM 

We consider a multibrane world scenario. The objective is to provide a model which 
concurrently solves both traditional fine tunings, the hierarchy and the cosmological con- 
stant problems. We first give the general multibrane solution and the 4-D effective action 
formalism. 



A. Multibrane World Solution 

Let us first recall the set-up of the static multibrane world. Consider N parallel 3-branes, 
located at y — Li,i — 0, ...N — 1. The brane at y = Li has brane tension o"j. Unless stated 
otherwise, the brane tensions are taken to be positive, except for branes at orbifold fixed 
points, for which we allow the possibility of negative brane tensions. In the uncompactified 
case, there are N + 1 AdS bulk spaces, with 5-D energy-momentum tensor T ab = Aig^ in 
the bulk between the (i — l)th and the ith branes. We start with the metric 

ds 2 = dy 2 + A(y) [-dt 2 + exp(2#t)<%dxW]. (36) 

Note that a rescaling of (t, x l ) rescales both A(y) and H, but leaves H 2 /A(y) invariant, so 
we have the freedom to fix the overall normalization of A(y). The Cos component of the 
Einstein's equation G a b = K 2 T a b = ^G^T a b (see for example Eq.(|2T|)) is satisfied trivially, 
while the Goo and the G55 components give, respectively, 
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A" 
~A 

2 



-) 
A 



2H 2 2k' 
— + - 

4ff 2 2k 2 



TV N-l 

Y, A i®(y - Li-i)®(Li -y)-J2 a ^(y - L i 

i=0 i=0 
N 

Y,^®{y-Li-x)Q{Li-y) 

i=0 



(37) 



(where — L_i = Ln = oo in the uncompactified case). The Gij components do not yield any 

additional equations. It is convenient to define = K 2 <Ti/3 and k{ = \J /t 2 Aj/6, where both 
have mass dimension 1. 

Integrating the Goo component of Einstein's equation across the % th brane yields the 
Israel jump condition ||23|| : 



lim 

y-*Li-\ 



~A 



— lim 



A' 
T 



= ~2qi. 



(38) 



In general, there is also an equation of motion for each brane; these are discussed for the 
5-form field strength case in Appendix For stationary branes only the 5th component of 
the equation of motion for the embedding coordinates is non-trivial. Since the derivatives 
of the metric at the brane are discontinuous, we have to average over the two sides of the 
brane. In the above metric ansatz, the equation for the ith brane reduces to 



lim 

y-*Li-] 



A'\ .. (A! 
— + lim — 

A J y^Li-\A 



(39) 



(see Appendix [A] for details). In the 5-form field strength model with M 5-form field 
strengths, Eq. (E2f) gives 



(4) + 



A iH 



A, 



(40) 



This means that, in the stationary situation, the average of A' /A of each brane is simply the 
difference of the bulk cosmological constants on its two sides. In the symmetric case, i.e., 
Aj + i = Aj, the limiting value of A' /A from the two sides are equal in magnitude but have 
opposite signs, as expected. Eq. ([39]) suggests that equations analogous to Eq. (flOD ought 
to hold for other models in which the bulk cosmological constant emerges as an integration 
constant. This is indeed the case for the generalized model fl3~4|), where 



' v f A '\ f A ' 
lim — + lim — 

y-*Li+\Aj y^Li-\A 



j 



+i 



A; 



(41) 



This completes the set-up of the problem. Not surprisingly, solutions of the Einstein equa- 



tion, Eq. (37), for stationary branes automatically satisfy the brane equations of motion, 
Eq. (p9|), when the jump conditions Eq. (p8[) are imposed. The approach we shall take is to 
first solve the Einstein equation (|37[) in the bulk, then impose the jump condition (|38D , and 
finally show that the brane Eq. (|39| ) of motion (for each brane) is automatically satisfied. 

Similar systems have been studied in Ref [p4}| . Defining ki = y / K 2 Aj/6, the solution to 
Einstein's equation (|3~7D in the bulks is 
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Mv) 



H 2 smh 2 [ki(y - yi )] 



k 2 



<y < Li 



(42) 



(as before — L_i = = oo in the uncompactified case). Continuity of the metric at each 
of the branes imposes the constraints, 



sinh 2 [fcj(Lj - y t )] sinh 2 [A;j +1 (Lj - y i+1 )] 



k 2 



and the jump condition Eq. (|38|) gives 

ki cosh[ki(Li - yi)] k i+1 cosh[k i+1 (Li - y i+1 ) 



(43) 



(44) 



smh[ki(Li - yi)] smh[k i+1 (Li - y i+1 )] 

It is easy to check, using both the continuity and the jump equations, that the brane equation 
of motion fl39|) is automatically satisfied for each brane. The 4-D cosmological constant as 
seen by observers on the ith brane is 

H 2 k 2 



H (Li 



(45) 



A(Lt) S mh 2 [k l (L l -y l )Y 

where A(Li) is the value of the warp factor at the ith brane i.e. A(Li) = A(y = Li). Hence 
there is an inherent sign ambiguity in the definition of H{Lj): 

■ ki 



H(Li 



(46) 



smh[ki(Li - yi)] 

Since k{ > by definition and we assume H(Lj) > 0, the sign we must choose is determined 
by the sign of sinh[/cj(Lj — y^)], or, equivalently, the sign of A' at y — > L^ . Inverting this 
expression, we can express the separation between two consecutive branes, A/, = fcj(Lj — 
Lj_i), in terms of the Hubble constants on those branes: 

Ak = -sU sinh- 1 ( ki ) + s- sinh" 1 ( Jp—] (47) 



\H(Li^) 

where sf_ x is the sign of A' at y — > Lf_ v 
(a) (b) 



H(Li 



A(y) 




A(y) 




A(y) 




L , 



i-1 



i-1 



FIG. 2 Schematic diagram of the possible types of behavior of the 
metric coefficient A(y) in the bulk between the branes. 
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The metric coefficient A(y) can have three different types of behaviors in the bulk be- 
tween the branes, as illustrated in Figure 2. We see from Eq. (|42]) that A(y) may decrease 
monotonically, increase monotonically, or must vanish at an intermediate y, which corre- 
sponds to a coordinate singularity of the metric, Eq. (|3(J), a particle horizon. Such horizons 
are discussed further in Sec. |Vll| . 

Because the definition of A/j incorporates both the brane separation and ki, these equa- 
tions can be thought of as determining (L, — Z/j_i) from H(Li) and H(Li-i) for given ki, or 
as determining ki from H(Li) and H(L i+ i), with the brane separations governed by some 
higher energy dynamics. In terms of the bulk cosmological constants on either side of the 
ith brane, 

H 2( L j = [gf ~ (fei+i + k i) 2 M - (fci+i - h) 2 } ^ 



in agreement with Ref. M5\. Note that Eq. appears to give the expansion rate on the 
ith brane in terms of local quantities - the bulk cosmological constants on either side of the 
brane, and its tension - but in the context of a full multibrane spacetime, these quantities 
may depend on the positions of the other branes, and the locations of particle horizons, 
which are not necessarily nearby. The form of the metric in Eq. (^) assumes a constant 
H, but when matter density p on the brane is included, one can easily solve for H 2 on a 
particular brane without solving for the metric in all regions of spacetime P,p6 ,p7,p5| . The 



result is simple: one simply replaces q% —>■ qi + n 2 p/3. This result provides a determination 
of the 4-D Newton's constant Gn by requiring H 2 to have the standard form || 

H 2 ^87rG N (A eff + p + ...)/3. (49) 

As we shall see later, this way of determining Gn is valid only when A(y) is peaked at the 
particular brane where the expansion rate is evaluated. 



B. The 4-D Effective Action 

Consider the 5-D action S^ 5 ' for the 5-D gravity plus scalar fields confined on branes: 



S^= [dyd 4 xV=g ^ + A + J2 S (y-L i )^(-a i + L i ) + / ^%^d 4 x (50) 

where the bulk cosmological constant A (A > for AdS) is piecewise constant, <7j is the 
ith brane tension, the surface term is the Hartle-Hawking term, needed for theories with 
boundaries, and 

U = \gr<kKh, - \m 2 Hi & - \4t -^2+"- ( 51 ) 

where the scalar field <pi is confined to the brane at L 4 . Recall that 87rG5 = k 2 = 87r/M 3 . 
Generically, the scalar field masses Mm (which are real after spontaneous symmetry break- 
ing) and the cut-off masses M ci are expected to be comparable to M in order of magnitude. 
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The characteristic scales brane tensions and the bulk cosmological constant are also expected 
to be set by M, i.e., Oi ~ M 4 while A ~ M 5 . The couplings A, are of order unity. 

Following Ref ||, we can integrate out the y direction to obtain an effective 4-D theory. 
As we shall discuss in more detail, the zeros of the metric A(y) are particle horizons. That 
is, for an observer on a brane, it takes infinite time for a light-like signal to travel from 
the brane to the particle horizon. This implies that the integration over y is only over 
the region between particle horizons that includes the position L of the visible brane. We 
can decompose the \iv components of the 5-D Ricci tensor R^J into the 4-D Ricci tensor 
R$(ffin>) an d the extrinsic curvature K^ v (and K = g^K^): 

R$ = R$ - 9 5 %K^ - g 55 K^K + 2g 55 K*K Xv . (52) 

where we recall the 5-D metric 

ds 2 = A (y) y^ (x) dx"dx u + dy 2 (53) 

so that the pullback on the Lj brane y^ v = g^ v = A(y = LAy^x), and \f^g = A 2 (y) \f—y. 
Using the fact that K^ u = g fMU A'/2A (where the prime indicates derivative with respect to 
y) and K = 2 A' /A, this gives 

J=jjR<® = ^ [A(y)RM - g^dyK^ — K 2 + 2K x ?K Xp + R 55 = 

^pf [A{y)RP> - {A{y)') 2 - AA(y)"A(y)] . (54) 

We substitute A(y) and its derivatives into the 5-D action S*- 5 -* and integrate over y to 
obtain the 4-D low energy effective theory. This "integrating out" of the 5th dimension in 
g(5) yields the effective 4-D action S*- 4 -*: 



= J d A xpR^/2K 2 N - J d 4 xpK f f + 
£ / d'xA 2 (LA yfj 



(55) 



2 A (L,) 1 T ' im 2 HY% m M 2 cl 
We see that the 4-D gravitational coupling is given by 

\h = h I A(y)dy - (56) 

Suppose that A(y) is peaked at the zth brane, namely the Planck brane, it is convenient to 
fix the normalization A(LA = 1; then, to a good approximation, we have 

G " - G *T§f- (57) 

">i ~r n/i+i 

for the region between particle horizons that includes the ith brane, and 87rG5 = k 2 . Col- 
lecting all the contributions to the effective cosmological constant A e jj, we have 
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A eff = J [-AA(y) 2 + {{A{y)'f + AA(y)"A(y))/2K 2 }dy 
-2A(y)' A(y) / K 2 \ boundaries + Y,A 2 (Li)a { 

i 

= - j AA{yfdy-^- 2 j{A'{y)fdy + Y, A2 ( L ^ ( 58 ) 

i 

where A" is singular at the branes. An integration by parts on the A" A term removes the 
Hartle-Hawking boundary term, as expected. The integration over y is between particle 
horizons and A is piecewise constant. Besides the contributions from the brane tensions and 
the bulk A, we see that there is a contribution from the 5-D Ricci scalar w^> to A e //. We 
can now re-express hi (or A) and yi in A e ff in terms of the brane tensions qi (or 04) and the 
brane separations Li. A priori, we expect K e ff to be of order Planck scale, but for large 
brane separations, we shall see that the boundary conditions at the branes fix the piecewise 
constant bulk A so that an almost exact cancellation among the terms in Eq. (|53) renders 
K e ff exponentially small. 



We now redefine the field fa = J A (Li) fa to absorb the factor A (Li) in the kinetic term 
of the Lagrangian Lj, so the effective 4-D low energy action becomes 



2 ^ e ff + * J CFT+ 




7 mfjiti u fa 

i,ti<Pi,v o \<Pi 2 + 

2 ml- 



(59) 



where m 2 Hi = A(Li)M^ and m 2 ci = A(Li)M 2 i. It is important to note that the Newton's 
constant and the effective cosmological constant A e ff are the same everywhere, so observers 
on the Planck (hidden) brane see the same Gn and A e ff as observers on the TeV (visible) 
brane, even if the visible brane tension is negative, as is the case in the Randall-Sundrum 
model. To be more realistic, L visib i e may be replaced by the standard model Lagrangian 
density. In S^ 4 \ we have also included a conformal field theory term he ft- Using the 
AdS/CFT correspondence |J, the effect of the Kaluza-Klein (KK) gravity modes may be 
incorporated into a conformal field theory on the brane. This should be the case when 
the model covers the region between particle horizons, where the gravity KK modes have 
a continuous mass spectrum. In the orbifold model, where the gravity KK modes may be 
discrete, in which case Lcft should be replaced by another appropriate strongly interacting 
field theory. 

Note that the warp factor A(y) does not appear in the 4-D effective action. Although 
we have chosen the warp factor at the Planck brane to be one, we could have chosen the 
warp factor at the visible brane to be one instead. The physics depends only on the ratio of 
warp factors on the two branes, which is unaltered by this rescaling. However, if we choose 
A(L V isibie) — 1; then the electroweak scale should be taken to be the fundamental scale, 
while the Planck mass is a derived quantity. 
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IV. TWO BRANE WORLD 



To render the general analysis more transparent, consider the simplest case, a two brane 
model. Setting L — 0, L x — L, the solution for the bulks outside the branes is 

A(y) = H2 ^ 2 My-yp)] {y < 0) 

Kq 

,2 



m= ^mz£k (0<y<L) 

My) = ^^'[M»-f)l (y > L) , ( 60) 

/Co 



where yo, Hi and ?/2 are constants, and give the locations of the zeros of A(y) when they are 
included in the domain of the spacetime in the fifth dimension. 
Continuity of the metric at the branes implies that 

sinh 2 (k y ) _ sinh 2 (£;i?/i) 

P ~~ P 

sinh 2 [fci(L - y ± )] _ sinh 2 [£; 2 (L - y 2 )\ ( , 

P P ' ^ ' 




FIG. 3 The two brane model in the uncompactified case, with the 
go brane at y = and the q\ brane at y — L. In the S 1 /Z 2 orbifold 
version (boxed), the two branes sit at the two fixed points. The 
metric factor A{y) is shown schematically. 



The jump conditions at the two branes are 
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ko cosh (ko yo ) k\ cosh(/c 1 y 1 ) 



9i- (62) 



sinh(/c 2/o) sinh(feiyi) 
fci cosh[A;i(L - y x )] k 2 cosh[k 2 (L - y 2 )} 
sinh[A;i(L - yi)] smh[k 2 (L - y 2 )} 

The expansion rate seen by observers on the brane at y = is H(0) = H/ ^ A(0), where 

H /uq k x 

A(0) smh 2 (k Q y ) sinh 2 (k 1 y 1 ) 

_ [fc 2 - (^o + go) 2 P 2 - (** - V 21 



(63) 

o 



with k\ — (ko ± go) 2 > or < 0, in agreement with Ref P5[] , which uses a slightly different 
approach. Similarly, the expansion rate seen by observers on the brane at y = L is H(L) = 



H/y/A(L), where 

H 2 k\ k\ 



A(L) smh 2 [k 2 (L - y 2 )} sinh^fc^L - y x )\ 
[kl-(k 2 + qi ) 2 }[k 2 -(k 2 - qi ) 2 } 

with k 2 -(k 2 ±q 1 ) 2 > or < 0. We can rescale t so that A(0) = 1, and the Hubble constants 



on the two branes are, respectively, H(0) = H and H(L) = H/^A(L). Note that although 
Eqs. ( |63"D and fl6l]) appear to determine the expansion rates on the two branes completely in 
terms of local quantities (i.e., the local brane tensions, and bulk cosmological constants just 
outside each brane), the values of these quantities on/near the two branes are connected via 
k\ and y x . 

The 4-D Newton's constant can be determined by introducing a small matter density 
p to the visible brane, that is, q — > q + n 2 p/3, finding the Hubble constant H, as in Eq. ([49]) 



and Refs. |8|,|25|-p7| , and then requiring that H 2 = (8^(^/3) (Agfj + p); the result is 

4:7iG N q = K 2 a k [l + 2K 2 A e// (2a + k )/3q Q k ] (65) 

where a = q — k . Although L— dependent, the correction term is small if GjyA e ff/k 2 <^ 1, 
which is the case here, so, to a very good approximation, we have 

K 2 a k 

G N = — . 66 

4irq 

Positivity of Gn requires ao > 0; to be specific, let us consider 

< h < q - fc . (67) 



Since the expansion rates H and H(L) given in Eqs. (|63|) and fl64|) depend on fci, our goal is 
to express H and H(L) as functions of L and the parameters k , k 2 , q and q x . This requires 
an expression relating L and k x ; from Eqs. (|6T| ) and ( |6^) we find 
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fclL = sinh " 1 (^)) +sinh "(^Z)) 

= sinn 

\^[k*-(k Q + qo y}[ki-{k,- q ,y} 

+ sinh- 1 f - 2hqi ). (68) 

K ^[k!-(k 2 + qi m!-(k 2 - qi y]J 

Here, Eq. (|68|) is regarded as a relation that determines ki in terms of L, q , q x , k and k 2 . 
In the 5-form field strength model, ki, k and k 2 will adjust together as the background field 
strength is determined as a function of L, q and q\. In general, H(0) ^ H(L). Because of 
the condition (p7|), H — > as fci — > cto = <7o — &o from below. This means that the expansion 
rate as seen by observers on the visible brane becomes exponentially small for large L, 

H 2 « W Q e- 2 ^ L - c ^ (69) 

where ao > 0, and sinh(C) = k\/H(L). As k\ — > ao from below, -ff(L) approaches a L- 
independent constant given by Eq. (|6"ID and so does C. This implies that A e ff becomes 
exponentially small as L increases. 

Note that there is only one integration constant for the bulk cosmological constants. 
Consider the model of Eq. (|34]). Here h% = —k 2 (J)/G, k\ = —K 2 ((f) + /i )/6 and k\ = —n 2 ((f) + 



Po + A*i)/6, where the integration constant is negative and the brane charges /io and n\ 
are essentially arbitrary constants. 

In the symmetric case, where k$ = k 2 = k and qo = qi = q (and let a = q — k = «o), we 
have H(L) = H(0) = H, and y\ = —yo = y 2 = L/2. The constant C in Eq. ( |69"D becomes 
C = aL/2, so, including a matter density p (which is treated as a perturbation), 

H 2 w Aa 2 e~ aL + ^!^ p . ( 70 ) 

Thus, A e ff still decreases exponentially with L, but slower than in the nonsymmetric case. 



A. The S 1 /Z 2 Orbifold Model 

To this point, we have concentrated on spacetimes that are noncompact in y, but similar 
results can be derived for the compactified case. First, we may choose to identify ko = k 2 = 
k\ and derive ki and H in terms of L. Next, we can compactify the y direction to a circle 
S 1 of length 2L. Placing the branes at y = and y = L, we may identify the two sides of 5" 1 
to obtain a line segment. That is, we perform a Z 2 orbifold, with one brane sitting at each 
of the two fixed points (y = 0,L). This S' 1 /Z 2 orbifold model is particularly simple, since 
there is only one bulk space between the branes sitting at the two end points. This may be 
considered as an expanding (non-supersymmetric) version of the Horava-Witten model PS 



and is discussed in Ref 0]. The model has branes with tension a = 3qo/K 2 at y = and 
o'l = 3qi/n 2 at y = L, separated by AdS space with bulk cosmological constant —6k 2 /n 2 , 
which is treated as an integration constant. The solution is 

A(y) = ^smh 2 [k(y-y )] (71) 
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for L > y > 0. Because of the orbifold symmetry of the model, the jump conditions are 



2k /q = tanh % 

2k/ qi = tanh[A;(L — y )] 

at y = and L, respectively. Combining the jump conditions implies 



tanh kL — q\/2k 



(72) 

(73) 
1 as 



The expansion rates on the branes are H(0) = H/yA(0) and H(L) = H/yA(L), where 
H 2 (0) = k 2 [(q /2k) 2 - 1] and H 2 (L) = k 2 [( qi /2k) 2 - 1]. We see that L is related to k via 



2k 1 - (q 1 /2k)t&nhkL 1 

if q\/2k = ±1, then qo/2k = =Fl irrespective of kL, but for q\/2k ^ 1, qo/2k 
kL — ► oo. According to our viewpoint, Eq. ( j73|) determines given q , q\ and L. 



L 



Eliminating k, we find that 



2k 



In 



1 + 2k/q ) (1 - 2fc/ 9i ; 
[l-2k/q )(l+2k/ qi ) 



(go + qi? - 2g git 2 - 2q\t 2 + (g + gi )A 
(go + gi) 2 - 2g gi* 2 + (go + ?i)A 



H 2 (0) = ^ 



where t = tanh kL and A = J (q + qi) 2 — Aq qit 2 



(74) 



(75) 



H 




FIG. 4 Hubble constant H 2 on the brane at Lq = as a function of 
the dimensionless separation I = kL and the hidden brane tension 
qi, in units such that the observable brane tension q = 1. Note 
that H 2 = for <7i = — 1, the RS scenario. 
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For large values of qoL, this reduces to 



h 2 (o) * go(gi + go) e - goL 

(9i - 9o) 



(76) 



where /c ~ go(l — e~ qoL )/2. Thus, for g £ 3> 1 and \qi\ > qo > 0, the cosmological constant 
on the y = brane is positive, and exponentially small. Moreover, although Eq. (|76|) may 
appear singular as gi — > q in fact 



H 2 {0) 



-q L/2 



(77) 



in that case. 

The behavior of Eq. (|75|) is shown in Figure 4 where we have taken q = 1 in some units. 
For gx > q (i.e., q\ > 1), H 2 = H 2 (0) brane becomes exponentially small for large / = kL. 
As qi decreases toward qo = 1, H 2 still drops off exponentially, but at a slower rate. For 
still smaller g^, H 2 goes to a constant ~ 1 for large / implying a large expansion rate on the 
L = brane. However, H\ = H 2 (L) then becomes exponentially small, as the roles of the 
branes at y = and L reverse (and we may then identify the brane at y = L as the visible 
brane). 



RS 




FIG. 5 Hubble constant squared versus separation. More precisely, 
er l = e~ kL vs H 2 , which is proportional to the 4-D effective cosmo- 
logical constant A e ff. The dashed line corresponds to our model. 
The solid vertical line corresponds to the Randall-Sundrum model. 
The horizontal dotted line corresponds to the single brane model. 
/ = 00 at the intersection of the 3 lines. 

If qi = —qo, as considered in 0, then y = ±00, so H(0) = H(L) = (see Figure 4) and 
therefore |go| — 2& for any finite non-zero L. In this limit, 



23 



A(y) = smh 2 [k(y - y )}/ sinh 2 [% ] - e" 2 ^, (78) 

the same as in the RS model ||. In general, our two brane model is qualitatively different 
from || , since it involves two branes with tensions of different magnitudes and not necessarily 
opposite signs, and interprets the bulk cosmological constant as an integration constant 
derived from the brane tensions and separation, but the RS model can be obtained as an 
appropriate limit of the more general model presented here. 

Figure 5, in which we plot e~ l versus H 2 , illustrates the key features of Eq. ([751). The 
dashed line represents the two brane model described above, and as I increases, the A e ff 
follows the dashed line and approaches zero. For the RS model q + q^ = = H 2 for any 
separation I. (Of course, g + <?L = is a fine-tuning.) Thus, the corresponding curve for the 
RS model is the solid vertical line at H 2 = 0. Suppose we start with only one brane, which 
may be considered as the I = oo limit of the two brane model. Then we find that the visible 
brane (at y — 0) can have any value of A e ff. This happens because the problem has reduced 
to that of a single brane, i.e., the jump condition for the brane at y = L = oo is absent, so 
k is no longer determined. It follows that H 2 , and hence A e //, can have any value; this is 
shown by the dotted line along the x-axis. This illustrates that the difference is due to the 
interchange of limits, i.e.,, the difference between solving the Einstein equation and before 
taking I — > oo, versus starting with I = oo before solving the Einstein equation. For large, 
finite /, we see that one does not have to tune the brane tensions to obtain an exponentially 
small A e ff, as long as the brane separation is relatively large. 

For positive tension branes, there will always be a zero of A(y) between the branes. 
This "particle horizon" , to be discussed below, is a barrier through which communication is 
impossible in finite time. In light of this, it is more useful to express the Hubble constant 
on the visible brane not in terms of the tension of the brane at L and its location which can 
never be detected, but only with q , and yo, the detectable quantities. Then Eqs. (|75| , |76D 
are 



H 2 (0) = ^ 1 -tanh 2 (% ) ^ ^e" 90 ^ (79) 



2 „2 



4 
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where yo is a function of q^, Ql and L. We discuss the role of horizons in determining the 
expansion rate on a given brane more completely in § |V 11 (J| . 

For positive tension branes, we shall see that the presence of the interbrane horizons 
prevents us from explaining the mass hierarchy problem, although we can obtain an expo- 
nentially small cosmological constant on one of the two branes - which we can identify as 
the visible brane - in that case. When one of the branes has negative tension, there are no 
horizons, and the warp factor will be small on that brane. Hence, we identify the negative 
tension brane as the visible brane, and can explain the mass hierarchy problem if the warp 
factor A(y) is sufficiently small on that brane. In that case, we also find that the expansion 
rate on the visible brane is small, but not small enough to explain the cosmological constant 
problem: when scaled appropriately, the expansion rate is of order the Higgs mass, which 
is far smaller than the Planck mass, but far larger than the expansion rate of our Universe. 
To find a solution to both problems, we need more than one separation distance and one 
bulk integration constant. In general, this may be achieved in a multibrane world, with 
more than two branes. However, we may also realize such a solution in the two brane case, 
provided that the brane world is compactified but not orbifolded. The properties of this 
model are discussed in detail in §[V], and from a different perspective. 
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V. THE COSMOLOGICAL CONSTANT AND THE HIERARCHY PROBLEMS 



The above two brane orbifold model is able to provide an explanation for either the 
hierarchy || or cosmological constant problem but not both. However, it is possible to 
solve both problems simultaneously in the multibrane world generally. In this section, we 
present details of a relatively simple model that solves both problems: a two brane model 
in which the fifth dimension is compactified, but not orbifolded. To capture the key feature, 
which is generic in a multibrane world, let us first present a toy model. 

A. Probing A One-Brane Orbifold Model 

It is easy to construct a variation of the above two brane S 1 /Z 2 orbifold model that 
foreshadows how the mass hierarchy and cosmological constant problems might be solved 
simultaneously in a two-brane model. Cut out the line segment jjq < y < L and remove 
the q% brane of the above S 1 /Z 2 model, so the resulting orbifold model has length y (with 
orbifold fixed points at y = and y = y Q ). The particle horizon is now at the orbifold fixed 
point at y — yo, and A(y ) = = A'(y ) at that point. There is only one brane (the q 
brane at the fixed point y — 0) in this S l /Z 2 orbifold model. H 2 is still given by Eq. (|79|) 
for q y ^> 1, which is exponentially small. Let us treat the go brane as the Planck brane, 
with A(0) = 1, and introduce the visible brane as a probe brane (i.e., with negligible brane 
tension, as in Ref j3U|) at y = L v , somewhere between the two fixed points. The warp (or 
hierarchy scale) factor on the visible brane is then A(L V ) ~ e~ qoLv , whereas A e ff oc e~ qoyo . 
Since yo > L v , the inequality @ follows. In this simple two brane model, we see that, if the 
hierarchy scale factor is exponentially small, the 4-D effective cosmological constant must 
be exponentially smaller. We can of course choose to keep the q\ brane for a three brane 
model. 

This model is strictly correct only if the visible brane tension is exactly zero, which is 
a fine-tuning. A more careful treatment is necessary for an arbitrary visible brane tension, 
even if it is very small. 

B. Two Brane Compactified Model 

Consider a two brane compactified model, which has a S 1 compactified 5th dimension, 
with circumference L 2 — L , as shown in the following figure. The brane at L has tension 
o"o and the brane at Li has tension ai, both of which are positive. Recall that go = ^ 2 o"o/3 
and gi = k 2 o"i/3. These two branes are separated by L\ — Lq on one side and L 2 — L\ on 
the other side of the circle, where L is identified with L 2 . Without loss of generality, let 
cr > (Ti > (g > gi > 0). Note that we require two integration constants, namely, the two 
bulk cosmological constants A x and A 2 , or equivalently, k\ and k 2 . 

For a piecewise constant bulk cosmological constant, we can introduce two 5- form field 
strengths, or alternatively, start with the model (|33D with only one 5-form field strength. Be- 
cause the 5th dimension is compactified, the brane charges, which are constant parameters, 
must add to zero: eo = Here k\ = — k 2 (0— (e( ) +eo) 2 /2)/6 and k\ = — n 2 (4> — e 2 )/2)/6, 
where the integration constants are the constant background field strength e( ) > and the 
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constant < 0. This allows us to treat fa and k 2 as integration constants to be determined 
by the 5-D Einstein equation. 




FIG. 6 The two brane compactified model, in which the hierarchy 
problem and the cosmological constant problem may be simultane- 
ously solved. L is identified with L 2 , so the circle has circumference 
L 2 — Lq. The brane at Lq (Li) is the Planck (visible) brane. The 
metric factor A(y) is shown schematically. 



The solution A(y) in the bulks is given by Eq. ( [42"D and H(LA is given by Eq. 
continuity conditions of the metric A(y) at the branes are ( [43l) 



The 



sinh 2 [fci(Li - y t )\ sinh 2 [/c 2 (^i - JJ2)] 



k'2 



(80) 



sinh 2 ^^ - y 2 )\ sinh 2 [£;i(Lo - yx)\ 



k'2 



and the jump conditions at the branes are 



q = k 2 coth[A;2(^2 - 2/2)] - fa coth[/ci(Lo - yxj\ 
qi = k,! coth[A;i(Li - y x )} - k 2 coth[fc 2 (L! - y 2 )}. 



(81) 



It is straightforward to obtain the large separation behavior. However, let us first get an 
overall picture of the model. The above conditions can be rewritten as 



k\ C oXk 2 \fa(L x -y x )\-\ 



L2 (coth 2 [fc 2 (L 2 -!/2)]-l)(l-t 2 ) 

K 2 7 



t 2 coth[k 2 (L 2 - y 2 )\y 



(82) 



k'2 



coth 2 [k 2 (L 2 -y 2 )\ - 1 =kf 



2(coth 2 [k 1 (L 1 -y 1 )}-l)(l-t 2 1 ) 



[l-h cothifaiL! - yi )\y 
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/ 4-i.n (t w\ i coth^^Li - yi)] - ti 
g = fc 2 coth[£; 2 (L 2 - 2/2JJ) - h: 



1 - ticoth[/ci(Li - yi)] 



7 u,rz \-|\ , coth fc 2 (L 2 -y 2 ) ~h 

qi = k x coth fei(Li - yi) ) - fc 2 - -TT, 77 yp 83 

1 - t 2 coth[/c 2 (L 2 - y 2 )J 

where ti i2 = tanh[/ci j2 A/!i i2 ], where AZi = ki(Li—L ) and AZ 2 = fc 2 (L 2 — Li). Let us consider 
the allowed regions in the (fci, k 2 ) plane. 




FIG. 7 The allowed regions (white) of "fc-space" for which there 
are consistent solutions with i7 2 (L 12 ) > and L 2 — L\ > L\ — L . 
The regions with small cosmological constant solutions are shaded 
dark. 



The allowed regions in Figure 7 are those for which the cotanh terms in Eq. (]81f) have the 
correct signs to correspond to the behavior of the scale factor A(y) in the bulk between the 
branes. Some of the details appear in Appendix |B[ Figure 7 shows the physically distinct 
allowed regions, where one demands L 2 — L\ > L\ — Lq without loss of generality. We see 
that there are two allowed regions: 

(1) the "one horizon" region for k 2 > k\ is bounded by the k\ = line, the qo = k± + k 2 
line and the Li — L = curve. In this region, 

k 2 (L 2 -L l )>k l (L 1 -L ) (84) 
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that is, the particle horizon is between L\ and L 2 segment of the circle, whose bulk 
cosmological constant has a larger magnitude (k 2 > ki). (The other k\ > k 2 > region 
has an identical description, and they are related by a simple interchange of the two 
segments of the circle.) 

(2) the "two horizon" region, where there is a particle horizon in each of the two segments 
of the circle. 

Let us consider the "one horizon" region where k 2 > k\ > 0. This region is bounded by 
the k\ = line, the q = k\ + k 2 line and the L x — L = curve. The L\ — L = curve 
is given by the equation k\ = k\ + qoQi, which can be derived from fl83l) . In this region, 
k 2 (L 2 — L\) > ki(L\ — L ), that is, the particle horizon is between the L\ and L 2 segment of 
the circle, that is, the segment with a larger magnitude bulk cosmological constant (k 2 > ki). 
This region is further divided into two subregions, where [L 2 — L{) < (L\ — L ) (the upper 
part) and (L 2 — Li) > (Li — L ) (the lower part). Because of the relationship between AZj 
and the behavior of the scale factors in the bulk, for Figure 6 we have 

Ah = hiLi-Lo) 

Al 2 = k 2 {L 2 - Li) 



=smh "(77^))- smlrl (-m 

, I ko \ , I kn 



= sinh 1 (777- 



C. Large Brane Separations 

It is easy to see that as H(Li) — > 0, Al 2 — > 00, but All can remain positive and relatively 
small if H(L ) = H(L 2 ) is smaller than H(Li). The brane at L\ will then be identified as 
the visible brane, on which A(L{) is large compared to H{L{) there, yet still exponentially 
small when compared to the fundamental (Planck) scale. 

Following the remarks above, we expect Al 2 to be large (t 2 — > 1). The boundary condition 
at L\ can then be satisfied if k\ — > 0, or coth[fci(Li — yi)] — > —1. The later solution 
will permit the two different scales required. The second boundary condition will require 
coth[fc 2 (-^2 — 2/2)] — > 1- Writing 

t 2 ~ 1 - 2exp[-2A/ 2 ] 
coth[A;i(Li - yi)} ~ -1 - rjx 
coth[k 2 (L 2 - y 2 )} ~ 1 + r] 2 , 

where i]i t2 are expected to be 0(exp[— Al 2 ]). To zeroth order in ^i j2 , the jump conditions 
are satisfied if 

h = ^(qo- qi), k 2 = ^(q + q 1 ). (85) 

We imposed go > qi which give ki > as required for consistency. The boundary conditions, 
to first order in r\^ 2 lead to 
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2 



Al 2 



We now demand that Alx is large, but it must remain smaller than AZ 2 (see Appendix [B]) 
so tx — 1 — 2exp[— 2A/i], which gives 

tf 2 (L ) = # 2 (L 2 ) ~ (g + qxfe-^ +Ah \ (86) 

Figure 6 shows schematically that the scale factor on the L\ brane is exponentially smaller 
than on the Lq brane: 



A(L ) H\Lx)~ e 



to lowest order. For a compactified 5th dimension, with both Alx ~ (go — Qi)(Lx — Lq)/2 
and Al 2 ~ (go + qi)(L 2 — L\)/2 large, we find that A(L\) and H 2 (L ) and H 2 (Lx) are all 
exponentially small. The results are 

A(Li) ~ e -(9o-gi)(ii-io) 

H 2 (L ) ~ (g + g 1 ) 2 e -[(90+9i)(i2-Li)+(go-9i)(Li-io)]/2 

H 2 {Lx) ~ (g + gi ) V [9o+9l)(L2 - Ll) - (9 °- 9l)(Ll - Lo)]/2 . (88) 
Since Li > Lq and go > gi, the warp factor A(Li) is exponentially small. 



D. Gn and A e j^ 

Let us consider further the situation as shown in Figure 6 in terms of the 4-D effective 
action. The y = L = brane is referred to as the Planck (hidden) brane and the y = Lx 
brane is referred to as the TeV (visible) brane. (Note that the visible brane tension may be 
of order Planck scale.) We adopt A(0) = 1, so A(Lx) = e~ 2Ah , where Alx ~ (go — qx)Lx/2. 
We shall take A(Lx) — 10 -32 , so the mass m H i = J A(Lx)M H1 of the Higgs field 0i on the 



visible brane is around TeV, while the mass mm = y A(0)M HO of the Higgs field O = 0o 
on the hidden brane is comparable to the Planck mass. In this model, the y integration is 
over the entire S l . Up to exponentially small corrections, we have 

f f H 2 1 1 

J A{y)dy = J — sinh 2 [k (y - y )]dy ~ — + — (89) 

with k = kx-i k 2 for each side of the y = L brane, so, up to exponentially small corrections, 
this integral over y yields the same Gn as that in fl57|). In terms of brane tensions, 



2q 



2 ,2 W 2- (90) 
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Using the result Eq. (|8q) for H , we find that the effective 4-D cosmological constant is 

A e// * ^0 + ^) e -A h -A h (Q1) 
(cr - cti) 

where All ~ (g — — L )/2 and AZ 2 — (qo + qi)(L 2 — Li)/2. In terms of fa and mm, 

Gn and Ae/y are the same for observers on both the visible (i.e., the TeV) brane at y = L\ 
and on the Planck brane at y = Lq = 0. Since if 2 = /CjyA e ///3, we expect the same Hubble 
constant for the two branes. However, Eq. (|86|) gives the ratio H 2 (L )/ H 2 (Li) = A(L\). 
This difference arises because Eq. fl36| ) is calculated in the (7^, 0«, Mffj) frame, while 
properties calculated from the effective action are in the (7^, fa, mm) frame. That is, 
in the (7^, fa, mm) frame, the Hubble constant seen by observers on the visible brane is 
simply H 2 = H 2 (L ). (More details on such A(Lj) rescaling later.) 

As pointed out in Ref. ||, the Higgs mass m Higgs as seen by observers on the visible 
brane is given by m 2 Higgs = A(L\)M^j. Using Eqs. ( p5| , p^ , p!25| ) , we see that the observed 
Higgs mass is exponentially small compared to the Planck mass for large brane separation 
L 1 = L 1 -L > 0, 

m m 99 s „ 2k 1 k 2 M 2 H A(L 1 ) = ±<kM 2 h ((jI - a 2 ) 87r((T0 _ (Tl)Ll/3M 3 

mhanck (h + h 2 )M* SM^o ' 1 1 

for o"o > <J\ > 0. Thus, we obtain a Randall-Sundrum like explanation of the hierarchy 
problem: an appropriate choice of the the warp factor can explain the hierarchy problem. 
At the same time <fi\ and other higher-dimensional operators become more important in the 
effective theory as their couplings are exponentially enhanced, that is, m c \ becomes TeV 
scale. 



E. An Inequality 

It is clear from the above that the mass hierarchy problem can be explained by this 
model if A(L\) ~ I0~ 32 . What is not so obvious, but shall be explained in § |V111| , is that the 
appropriate statement of the cosmological constant problem is that H 2 (L )Gn ~ 1CT 122 . 
This may be understood most easily in terms of the effective theory introduced in §111 B 



where we show that, with appropriate choice of units, Gn and H may be regarded as 
the same on all branes, but the Higgs mass on any particular brane is related to Gn by 
GN m mggs ~ A(Li). Assuming that qo,qi ~ G N 1 ^ 2 , we see that Eq. ( p8|) (as well Eqs. ( |87D 
and (0)) implies that if we can explain the mass hierarchy problem, then the observed 
expansion rate, in Planck units, is smaller than the observed Higgs mass in Planck units by 
additional powers of A(Lj) 1, because A/ 2 > A^. Although the static models do not fix 
the size of expansion rate precisely, even when the ratio of Higgs and Planck masses is taken 
from observation, they guarantee that the value is far smaller than the Higgs mass. These 
considerations are based on the one- horizon model, and the inequality Al 2 > Al\. 

In the above analysis, we show how two exponentially small factors can be generated, 
one for the hierarchy problem and the other for the cosmological constant problem. As an 
input, we choose a large radius (L 2 — L )/27r for the compactification. It is important to 
find a dynamical reason for such a large radius. Here let us address a much more modest 
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question: why should the exponent for the effective cosmo logical constant, m Planck / A e //, be 
bigger than the exponent for the hierarchy factor, m 2 Planck j 'm 2 Higgs . Since both exponential 
factors have similar origins, we actually have an inequality between these two exponents 
that is satisfied in nature. 

The solution we have considered is the "one horizon" region in Figure 7. As we approach 
go = ki + k 2 and q± = k 2 — hi, both L\ — L and L 2 — L\ become large. Since Al 2 > Al\ 
(that is, (go + qi)(L 2 — L\) > (g — gi)(-^i — -^o)), we have the inequality (|J), namely, 
\n(m Planck / 'A e f f) > \n(m 2 Planck / 'm 2 Higgs ) , as is the case in nature. However, we can actually 
make a stronger inequality, that is, L 2 — L 1 > L x — L with g + q x > q — qi {i.e., the 
particle horizon is in the larger segment of the circle). Let us assume that the positions of 
the branes are fixed, that is, the stationary branes satisfy their corresponding equations of 
motion. There are two distinct cases here: 

(a) L 2 — L\ > L\—Lq > and k 2 > k\\ that is, the particle horizon is in the larger segment 
of the circle between the branes. This case is schematically shown in Figure 6. 

(b) The particle horizon is in the shorter segment of the circle. There are two equivalent 
descriptions of this case. We can either fix L 2 — L\ > L\ — L$ > 0, which implies that 
Al 2 < Ali and k\ > k 2 ; or choose L\ — L > L 2 — L\ > in the k 2 > k\ region with 
Al 2 > Ali. The Hubble constant in this second case (with the particle horizon in the 
shorter segment of the circle) is given by 

H\L ) = (q + qi ) 2 e-^ + ^ (93) 

where Al 2 = (q - q{)(L 2 - Li)/2 and AZi = (q + qi)(L l - L )/2. 

Nature will pick the lower energy (or least action) solution, that is, the one with smaller 
A e ff. Since the Newton's constant is essentially the same in both cases more discussion on 
this later), we can simply compare this Hubble constant to that in Eq. (j36|): 

(?o + qi)(L 2 - Li) + (g - gi)(£i - L Q ) > (g - qi)(L 2 - L t ) + (g + q\){L x - L ) 

=> L 2 -L x > Lx - Lq. (94) 

This implies that the particle horizon should be in the larger segment between the branes, 
as schematically shown in Figure 6, where L 2 — L\ > L\ — Lq > and k 2 > k\. Note that 
this inequality L 2 — L\> L\ — L together with g > gi > is stronger than the inequality 
(<?o + qi)(L 2 — Li) > (g — qi)(Li — L ). Using this inequality (|94]) and g > gi > 0, we can 
now compare the exponent in Eq. (BH), where 



2 

m ^ i99S ~ A(Lx) ~ e -(*o-<n)(Lx-L ) (95) 

m Planck 



and the exponent in Eq. (^6|), where we can write in terms of A 



effi 



A 2ffo(g + gi) e [( go+gl )(£ fl -xi)+( g0 - ft )(L 1 -x o )]/a / g6 ^ 

{qo - gi) 



to obtain 
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j > ln ( Kla^cA (g7) 

V A e// J \m 2 Higgs J 



as is the case in nature, where \ia(m% lanck / A e ff) ~ 2.3 x 122 and ^{mp lanck /m 2 Higgs ) ~ 
2.3 x 32. To get a feeling for the magnitudes of the various quantities, we see that choosing 
L 2 — Li ~ 2(Li — L ) and cr ~ 2a\ gives the correct order of magnitude. If we assume that 
the brane tensions are of Planck scale, the values of the cosmological constant and mass 
hierarchy observed in nature follow if (Li — L ) ~ ^ m ~p} a nck- Smaller brane tensions imply 
larger {L x - L ). 

In this model, the branes are stationary (that is, they satisfy the brane equations of 
motion) and both L\ and L 2 are stable. But we still need a dynamical explanation for 
why Li (or L 2 ) is large. However, if we understand the hierarchy problem via some means, 
then this brane world model provides an explanation of the cosmological constant problem. 
Now, the hierarchy problem may be understood from the renormalization group flow in 
4-D quantum field theory. It is well-known that the running of the couplings (marginal 
operators) are logarithmic versus the running of the mass terms (relevant operators). A 
well-known example is the unification of the couplings in some supersymmetric versions 
of the standard model. If we accept some version of this picture as a resolution of the 
hierarchy problem, then the cosmological constant problem is explained by the inequality 
(|j) that appears in this two brane world model. In fact, the running of the couplings and 
masses in 4-D quantum field theory can be recast as the holographic renormalization group 
flow in the brane world scenario |31|| . So, in a brane world model where the hierarchy problem 
is solved, the inequality (|j) tells us that we have an exponentially small 4-D cosmological 
constant. 



F. Two-Horizon Region 



There is another possible solution, namely, the "two horizon" region shown in Figure 7. 
In this case, the two branes are separated by particle horizons. Since we are approaching 
Qi = k 2 + k±, the value of H 2 /A(0) is finite, so the q\ brane again has exponentially small 
cosmological constant. Within the 5-D classical Einstein theory, the two solutions are topo- 
logically distinct and so we cannot compare them. However, one may imagine a deviation 
from the pure AdS situation, or in quantum gravity or string/M theory, where they may be 
compared. Which solution does the two brane system pick? As is discussed in more detail 
in Appendix [B[ the AD effective theory has a well-defined A e //, so we may compare the 
values of A e ff in the two cases. For large brane separation L 2 — L\, we have seen that A e ff 
is exponentially small in both cases. By comparing A e ff in the two cases, the "one horizon" 
case (where the hierarchy problem may be solved) has a smaller A e ff, due to smaller k 2 and 
maybe also k\. For finite brane separations, and with appropriate choices of qi, it is possible 
that the second case is preferred. One can imagine that the brane world starts with (two) 
particle horizons separating the two branes. As the brane separations (and the radius of S 1 ) 
increase, the bulk cosmological constants on two sides jump to the values for the hierarchy 
solution phase. We do not expect this to happen in the 5-D theory. However, within the 
string/M theory framework, such a transition is akin to a topological change, a possibility 
that deserves further study. 
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VI. MULTIBRANE WORLD 



It should be clear that a three brane or multibrane model can provide a similar solution to 
both fine tuning problems by the same mechanism; on the visible brane H 2 ~ exp(— 2A/ 2 ), 
with a coefficient dependent on all three brane tensions, and A(Li)/A(L Q ) ~ exp(— 2A/!). 
However, the 5-form field strength mechanism to provide piecewise constant bulk energy 
densities fails in an orbifolded space and it is not currently known how to achieve the 
different adjustable bulk cosmo logical constants. A multibrane world with 5-form field 
strengths providing the piecewise bulk cosmological constant can give a solution to both 
fine-tuning problems, in non-orbifolded spaces. Other variations of the models studied here 
should be explored further. 

We can now give a general picture of the physics when the brane separations are large. In 
the stationary situation, we shall treat the brane positions Li as fixed, while the parameters 
ki between branes and i/i are to be determined by the jump conditions, via Eqs. ([43|) and 
By solving the jump conditions, the Hubble constant seen by observers on the ith 
brane can be expressed in terms of the 5-D gravitational coupling k, its brane tension <jj 
and the tensions <7j_i and o"j + i and positions and Lj + i of its neighboring branes. For 
large brane separations the relation ([P]) between the bulk parameters ki and ki + \ and the 
brane tension becomes 

qi = ki coth[ki(Li - y { )] - k i+1 coth^+^L; - y i+1 )], 

-> s^ki - s?k i+1 (98) 

where sf is the sign of A' at y — > Lf, and the corrections are exponentially small for large 
separations. Therefore, for any of the four possible sign choices we see that H 2 (Li) becomes 
exponentially small. The implications of the various sign choices for the cosmological con- 
stant and hierarchy problems can be summarized qualitatively in the following table, where 
it is assumed that k i+ i > ki > 0. (The opposite inequality k i+ i < ki does not yield any new 
possibilities, as they can be obtained from the table below by exchanging some of the cases.) 
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Qi < -{k+i + ki) 




? 


+ 


+ 


V 


4 
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qi > -(h+i ~ ki) 




X 


+ 


+ 


V 



We consider cases for which H 2 (Li) > only. For case (1) the scale factor is locally 
peaked at the brane. The brane in case (3) can have no zeros between it and either neigh- 
boring brane, while the remaining cases may have a horizon on one side only. The qi column 
gives the sign of the brane tension; stability is expected for a brane with a positive tension. 
Generically, a negative tension brane is unstable, as in case (4), but, it is possible that some 
stabilization mechanism may be found in a more complicated model. Stability is assumed for 
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a negative tension brane if it sits at an orbifold fixed point, so that the fluctuation modes 
that will cause destabilization have presumably been projected out. In this case, the Z 2 
symmetry requires ki + % = k{. This is represented by " ? " in case (3) in the table. 

As we shall explain later, the Newton constant is always positive and is the same on all 
branes. The value of the 4-D Newton's constant Gn may be determined by the introduction 
of a small matter density in case (1), where A(y) is peaked at the brane, as given in (fl9|). 
However, the same method cannot be applied to branes in the other cases, where A(y) is 
not peaked at the brane (see Appendix |D|). 

In the above framework, and without fine-tuning, we see that there can be models that 
simultaneously 

(1) have a stable positive tension brane, 

(2) have an exponentially small 4-D positive cosmological constant, and 

(3) yield a solution to the hierarchy problem. 

Self-consistent models with all of these desirable properties require a combination of case 
(1) and case (2), so the simplest model must contain at least two branes. Generically, 
a multibrane model can accommodate both an exponentially small cosmological constant 
and a solution to the hierarchy problem. The model of Ref |3J corresponds to the case 
(1) in the Table, in which the hierarchy problem is not solved. We shall show that the 
hierarchy problem in this model is simply the standard one, so it may be solved by more 
conventional means. The Randall-Sundrum solution to the hierarchy problem corresponds 
to the converging limit of case(3), where ki + \ = hi, with a negative brane tension. The 
hierarchy problem may be solved when the metric is non-vanishing before it reaches the 
next brane. This is indicated by a " a/ " in case(2) in the Table. The simplest realization 
of such a solution to both the cosmological constant problem and the hierarchy problem is 
a two brane model in the compactified (but not orbifold) case. In this model, the hidden 
Planck brane corresponds to case(l) and the visible brane corresponds to case(2). 

In fact, the two brane compactified model, although simplest, is rather specific, and we 
can generalize the main result to prove that H 2 / m 2 Higgs <C 1 in a multibrane model. Suppose 
that the visible brane is the ith brane, with a warp factor A(Lj) relative to the value on the 
Planck brane, which is somewhere to its left, with no particle horizons in the intervening 
region. Let D{ = Li — Lj_i be the distance to brane i — 1 just to the left of the zth brane, 
where the warp factor is A(Lj_i) > A(Li), and let A+i = L i+ i — Li be the distance to the 
brane i + 1 just to the right of the visible brane, where the warp factor is A(Li + i) < A{Lj). 
(This presumes that there is no particle horizon between branes i and i + if there is, then 
take A(L i+ i) = and let D i+ i be the distance from the ith brane to the horizon.) 

Then instead of Eq. (B3), we find 



H 2 « Aq 2 A(Li 



A(Li 



A(Li_ 



-2qiD i+1 



InjAjL^/AjL^ 2 



x 



_ [A(L i+1 )]V 2 [A(L^)] D ^ g l . 

[A(L0] 1/2+Di+l/2Di y 

once again, we see that H 2 /m 2 Higgs is exponentially small, although the model is more compli- 
cated than before. (Consistency requires non-negative H 2 .) The two brane, uncompactified 
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version of Eq. (|99|), which is appropriate when horizons surround two branes, one of which 
is the visible brane (with tension qy and warp factor Ay) and the other the Planck brane 
(with tension qp and warp factor Ap = 1) is 

(100) 

where D is the distance from the visible brane to the horizon to its right, and Dp is the 
distance from the visible brane to the Planck brane to its left. Eq. ( j88|) can be obtained from 
the uncompactified model by assuming a periodic sequence of identical two brane regions 
bounded by horizons. 



H 2 



±q v A l J D « IDp e- 2 ^ D « 



2q v Dp 



VII. GLOBAL STRUCTURE OF THE SOLUTIONS 

All of the solutions discussed in this paper consist of regions of 5-D anti-deSitter space 
joined across the 3-branes, together with some additional identifications in some of the 
solutions. However, the form fl36|) of the 5-D anti-deSitter metric used here is unconventional, 
and in addition in some cases possesses coordinate singularities where the metric coefficient 
A(y) vanishes. Coordinate singularities of this type are present in all solutions for which 
all the branes have positive tension fj As discussed in the introduction, these zeros of A(y) 
are particle horizons, in the sense that signals traveling from a brane to the horizon are 
perceived by observers on that brane to take an infinite time to reach the horizon (although 
the elapsed proper time or elapsed affine parameter for null signals is finite). These horizons 
are closely analogous to the horizons in Rindler space. In this section, we show that the 
coordinate singularities are removable by exhibiting the coordinate transformation between 
the form (^) of the metric and a more standard form of 5-D AdS space. We also derive the 
global structure of the various bulk solutions we have been discussing. The key result that 
we shall obtain is that the solutions typically split into different disconnected components 
across the particle horizons. The global structure of the RS model has been discussed in 
Ref. ra. 



A. Different coordinate systems on 5-D anti deSitter space 



We start by discussing the global structure of the spacetime ( |3~6|) when no branes are 
present, so that Lj = — oo and L i+1 = +oo. The metric flSTf ) can be written 

ds 2 = dy 2 - A(y)dt 2 + A(y)e 2Ht [(dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ] , (101) 

where the coefficient A(y) is given by Eq. 



3 In the two brane case, the coordinate singularities can be avoided (by choosing the free param- 
eters appropriately) when the two branes are not identical. However, in the symmetric two brane 
case, a coordinate singularity between the two branes is unavoidable. 
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H 2 

A (y) = -p- sirm2 [h(y - yd] ■ 



We define new coordinate T, Y and X 1 for 1 < I < 3 by 

T = i tanh fofo - W )/2] {(1 + tf 2 x 2 )e m - e"**} 

F = i tanh [k^y - 2/i )/2] {(1 - H 2 * 2 )e m + e~ m } 
X< = #tanh [h(y- yi )/2] e m x\ 



where x 2 



x 



1\2 



X 



2\2 



+ (x 3 ) 2 . The inverse of this transformation is given by 



tanh 2 [ki(y - yi )/2] 
exp [2Ht] 

Hx l 



X 2 + Y 2 - T 2 

(T + Y) 2 
X 2 + Y 2 - T 2 
X 1 



T + Y' 



(102) 

(103) 

(104) 
(105) 

(106) 
(107) 

(108) 



where the sign of y — yi is taken to be the same as the sign of T + Y, and where X 2 = 
(X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 . In the new coordinates (T, Y, X 1 , X 2 , X 3 ) the metric takes the 
standard, conformally-flat AdS form 



4/Af 



;i + T 2 - X 2 - Y 2 f 



dT 2 + dY 2 + (dX L ) 2 + (dX 2 ) 2 + (dX 3 ) 2 



(109) 



See Appendix |C| for details [] 

We now discuss some properties of the coordinate transformation. From Eq. (|106|) , the 
original coordinates (t, y, x 1 , x 2 , x 3 ) cover only the region F 2 + X 2 > T 2 of the full spacetime. 
If we define R = VX 2 + Y 2 , the metric can also be written as 

4A 2 



ds 1 



(1 + T 2 -R 2 ) 2 



-dT z + dR z + R z dQi 



111) 



where dfl 2 is the metric on the unit three sphere. The conformal factor diverges at the 
hypersurface R = \/T+ T 2 , which is the timelike boundary at infinity of the spacetime. The 
Penrose diagram is Fig. 20 (ii) of Ref. p3 . 



We note that another conformally flat coordinate system is obtained from the transformation 



-m 



-m 



Y 



H tanh[ki(y - yi)Y ' H sinh [ki(y - y,)] ' 

In the coordinates (T, Y ) the metric takes the form 



ds 2 



1 



BY 2 



-dT 2 + dY 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 



„2\2 



„3\2 



(110) 



The domain of the original coordinate system (t, y, x , x , x ) is the region |T| > \Y 
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FIG. 8 The shape of the hypersurface a = ao in the new coordi- 
nates (T, Y, X 1 , X 2 , X 3 ). Two of the coordinates, X 2 and X 3 , are 
suppressed. 

Another useful coordinate system is given by R = a cosh?/;, T = asinhip, in which the 
metric takes the form 



a 2 ) 2 



-a 2 dip 2 + da 2 + a 2 cosh 2 1\) d£l\ . (112) 



The range of these coordinates is < a < oo and — oo < ip < oo. This coordinate system 
covers only the region R < \T\ of the full spacetime. From Eq. ( |10(j| ), the coordinate a is 
related to the original coordinates t, y, x l by 

a= |tanh[^(y- m )/2]|. (113) 

Consider now the hypersurface a = ao, where ao > is a constant. In the coordinates 
(T, Y, X 1 , X 2 , X 3 ), this hypersurface is the hyperboloid 

Y 2 + X 2 - T 2 = a 2 . (114) 
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From Eq. ([L12 ), the induced metric on this hypersurface is 

-dtp 2 + cosh 2 t/j dn 2 3 



< 4 W 



4«o 



(115) 



which is the spatially compact (k = +1), geodesically complete version of 4-D deSitter space. 
Using Eq. ( |113Q , we see that in the original coordinates (t, y, x 1 , x 2 , x 3 ), the hypersurface 
( |114|) consists of the union of the two surfaces 



and 



y = yi + 2tanh 1 (a Q )/k i 



y = yi - 2tanh 1 (a Q )/k i , 



(116) 



117) 



on the right and left hand sides of the coordinate singularity at y — y^. Thus, the two 
hypersurfaces (|116|) and ( |117|) are in fact connected, even though they appear disconnected 
in the original coordinates (t, y, x 1 , x 2 , x 3 ). The hypersurfaces ( |116| ) and ( |117|) are mapped 
onto the regions Y + T > and Y + T < respectively of of the hypersurface (|114|) . 
The hyperboloid a = «o is illustrated in Fig. 8. Note that the induced metric on each of 
the hypersurfaces (|116|) and (|117|) is, from Eq. (|3~6"D, the spatially non-compact {k = 0), 
geodesically incomplete version of 4-D deSitter metric [33[. Note also that it follows from 
Eq. ( |103| ) that if the vector d/dt is taken to be future directed for y > y^, then it is past 
directed for y < y^, just as in Rindler spacetime. 



B. The class of geodesically complete multibrane solutions 

The various multibrane solutions discussed in the earlier sections of this paper are 
geodesically incomplete. In particular, all the "branes" discussed in the earlier sections are 
surfaces of constant y, and it is clear from the discussion above that a surface of constant 
y comprises in reality just half of a brane. Our goal in this section is to find geodesically 
complete extensions of those multibrane solutions, in which the "other halves" of all the 
branes are present, and also to elucidate the global structure of those extensions. 

In attempting to find the geodesically complete extensions, it is convenient not to start 
with the multibrane solutions discussed earlier in the paper, but instead to start by directly 
constructing a class of geodesically complete multibrane solutions via cutting and pasting 
regions of 5-D AdS space in the form (|1 1 1|) . The metric is 



ds 2 = kHl + T2 _ R2)2 [-dT 2 + dR 2 + R 2 dn 2 3 \ , (118) 

where we have written k for fcj. Consider now a submanifold of this manifold given by 

a\ < R 2 - T 2 < a\. (119) 
We shall make use of three types of such submanifolds: 

• Type I, where a\ = — oo and < a\ < 1. Such regions contain the "particle horizon" 
R — \T\ and are spatially compact. They have a single 4-D boundary at a = a 2 - 
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• Type II, where < a\ < of < 1, which are spatially compact. They have two 4-D 
boundaries at a = a.\ and at a = a 2 . 

• Type III, where < a\ < 1 and a\ = 1, which are spatially non-compact. They have 
a single 4-D boundary at a = a%. 

Each such submanifold is characterized by three numbers: the values of a\ , a\ and of k. 

Consider now the class of solutions to Einsteins equations that one can obtain by gluing 
together submanifolds of the above type at their boundaries. It is clear from the number of 
boundaries that each type of submanifold possesses that such solutions, if connected, must 
be in one of the following two classes: 

• Class A: The solution starts at one "end" with either at type I region or a type III 
region, then has one or more type II regions (with either orientation, i.e, a can be 
increasing or decreasing), and then finishes at the other "end" again with either a 
type I region or a type III region. 

• Class B: The solution has no "ends" and instead consists of one or more type II regions 
joined together in a circular fashion. 

Additional solutions can of course be generated from these classes of solutions by performing 
additional identifications, and also by considering solutions with more than one connected 
component. 

The junction conditions in any of these solutions are the same as those obtained earlier 
in the paper, except that they are now expressed in terms of different quantities. Continuity 
of the induced metric ( |115[ ) across any boundary or brane requires that the quantity 



take the same value on both sides of the brane. The Israel junction condition is the require- 
ment that the rescaled brane tension q at any boundary is the sum of two terms, one from 
each submanifold on either side of the brane. The absolute value of each of these terms is 



The sign of each of these terms is positive for the boundary of a type I region and for the 
second boundary a = a 2 > cti of a type II region. It is negative for the boundary of a type 
III region, and for the first boundary a = cti < a 2 of a type II region. 

If one requires that all brane tensions be positive, then class A solutions with type III 
regions on both ends are disallowed, as are all class B solutions. Thus, there are two different 
types of solutions for positive tension branes: spatially compact class A solutions with type 
I regions on both ends, and spatially non-compact class A solutions with a type I region on 
one end and a type III region on the other end. 




(120) 




(121) 



2a 
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C. Relation to earlier solutions 



Consider now how the above geodesically complete multibrane solutions are related to the 
multibrane solutions discussed earlier in the paper, which were characterized by parameters 
hi, Li and y € . 

First consider the case of solutions with no particle horizons (which requires negative 
tension branes if the bulk regions are all AdS as in this paper). For these solutions, there 
is a simple and one-to-one correspondence between the solutions of Sec. [I l| and those of 
Sec. [VII B| above. For a given solution (ki,yi,Li), one can obtain the geodesically complete 
extension as follows. For each interbrane region L^_ x < y < Li where Lj_i and Li are both 
finite, one constructs a region of type II with the corresponding parameters obtained from 
Eq. ( |113j ), namely a± = min(&i,d2) and c*2 = max(di,&2), where 

a\ = tanh 2 [k^L^ - yi )/2] (122) 

and 

&l = tanh 2 [ki(Li - 2/0/2] . (123) 

For any regions Lj__i < y < Li where one of the boundaries is at \y\ = oo, one similarly 
constructs a region of type III. Then, one glues the regions of type II and III together to get a 
geodesically complete solution. The junction conditions will be automatically satisfied. Each 
interbrane region in the original solution maps onto the subset Y + T > of a type II or type 
III submanifold of the full solution. Each y = constant half-brane in the original solution 
acquires an additional half-brane in the full solution to make it a geodesically complete 
brane. 

Now consider solutions with particle horizons. First, for interbrane regions not containing 
particle horizons one follows the same procedure as above to construct the corresponding 
portions of the maximally extended solution. Second, for each interbrane region Lj_i < y < 
Li containing a particle horizon y = yi, there are two possibilities for the global structure. 

• One can map the entire interbrane region, including both sides of the particle horizon, 
into a single region of type I. Because the type I region has a single boundary a = ct2, 
this structure requires that the surfaces y = Lj_i and y = Li both map onto the 
boundary a = ct2, and therefore both surfaces must be equidistant from the particle 
horizon at y = yi. More generally, the entire multibrane solution must be reflection 
symmetric under the mapping y — yi — > —(y — y^, since the regions y > yi and y <y% 
are the Y + T > and Y + T < portions of the same bulk region. 

• If the above property of symmetry under reflection through the particle horizon is not 
satisfied by the multibrane solution (ki,Li,yi), then the only possible global structure 
is as follows. One identifies the region y > yi to the right of the particle horizon with 
the Y + T > portion of a type I submanifold, and one identifies the region y < yi 
to the left of the particle horizon with the Y + T < portion of a different type I 
submanifold. Therefore the maximally extended multibrane solution can have several 
disconnected components. 
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The general structure of a multibrane solution (ki,Li,yi) without any additional iden- 
tifications and without any special symmetries is therefore as follows. The spacetime can 
be split up into cells bounded by particle horizons, but with no particle horizons within 
any given cell. Then, the maximally extended solution consists of a number of different 
connected components, one for each cell. Since the physics of each connected component is 
independent from all the others, it is natural to restrict attention to solutions with only one 
connected component or cell. 



This global structure justifies our prescription of Sec. [Ill B| for obtaining an effective 
4-D theory of restricting the integral over the fifth dimension to between pairs of particle 
horizons. The global structure also explains the fact that the Hubble constant H(Li) on 
any brane can be expressed entirely in terms of the tensions and locations of the branes in 
the same connected component, without reference to the detailed arrangement of external 
branes. For example, consider two branes in an uncompactified spacetime. Let the two 
branes be sandwiched between two horizons, and assume there is no horizon between the 
branes. Let the two branes be at y = and L, and let the bounding horizons be at y = — yo 
and y = L + y 2 . Assume that A(0) — 1 > A(L). The expansion rates on the two branes 
turn out to be 



H 2 (0)^4q 2 L 
H\L) « Aql 



1 + 



2q L L 
ln(lA4(L)) 
2q L L 



l+V2/L e -2q L y 2 

[A{L)] y2/L e- 2qLV2 , (124) 



[A(L)] 1+y ^ L e 

2 



where qi is the tension on the brane at L, and we have written the solution in terms of 
A(L), which may be regarded as an observable quantity, the ratio {mmggsl 'mpianck) 2 found 
on the visible brane. 

For solutions with additional identifications, the existence of particle horizons can be 
compatible with the maximally extended solution having only one connected component. 
For example, consider the the two brane compactified model of Sec. |V[ That model has 
branes at y — Lq, y — Li and y = L 2 , but L and L 2 are identified. When there is only one 
horizon, then the maximally extended solution is connected and starts with a type I region 
at one end, then has a type II region, then has another type I region at the other end. When 
there are two horizons the maximally extended solution has two connected components. 

The global structure analysis also shows that the multibrane solutions of the earlier sec- 
tions are a subset of the full set of multibrane solutions, for the following reason. In the 
earlier sections, it was always assumed that the value of the bulk cosmological constant was 
the same on both sides of any particle horizon in a given interbrane region. However, since 
the two sides of a particle horizon typically correspond to two different type I regions in the 
maximally extended solution, there is no reason to make this restriction. The most general 
class of multibrane solutions, when translated into the old (t, y, x 1 , x 2 , x 3 ) coordinates, allows 
interbrane regions with particle horizons where the value of k appears to change discontin- 
uously as one crosses the particle horizon. The discontinuity is only apparent as the two 
sides of the particle horizon are not physically connected. 
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VIII. 4-D NEWTON'S CONSTANT AND PARTICLE HORIZONS 



There are some issues that we like to resolve in this section: the determination of the 4-D 
Newton's constant Gn and when a large warp factor helps to solve the hierarchy problems. 
We have already assumed the some of the results of this section in earlier discussions. Here 
we like to clarify and justify them. 

A. The Issue 

Naive determinations of Gn sometimes yield different answers. In the multibrane model, 
we face the following issue. We know of 3 approaches to determine the 4-D Newton's constant 
Gn'- 

(i) integrate over the 5th dimension in the 5-D action to obtain the low-energy effective 
4-D action ||, as discussed earlier; 

(ii) solve for the trapped gravity mode and determine its probability (wavefunction 
squared) on the brane 0; 

(iii) calculate the Hubble constant and determine the coefficient of the 4-D cosmological 
constant or the matter density contribution term ||. 

So far in this paper, we have used the approaches (i) and (iii). In approach (iii), note 
that the value of Gn depends only on the local (i.e., neighboring) properties of the brane 
on which we measure the value of Gn'- 

G * = G *irz¥- ( 125 ) 

for the ith brane, where 87t(j5 = k 2 . Recall that fcj (fcj+i) is a function of as well as g,_i 
and Li — L^i (qi+i and L i+1 — Lj) (where the explicit dependence is model-dependent), so 
Gn depends on neighboring properties only. Specifically, it does not depend on the total 
number of branes in the model. On the other hand, if we use the first approach, and naively 
integrate over the whole 5th dimension (i.e., from y = — oo to y = oo in the uncompactified 
case), we will get a Gn that depends on properties of all the branes in the model, specifically, 
on N, the number of branes. This is also the case in the second approach. If we normalize 
the gravity trapped mode over all y, we will typically get a Gn that is roughly N times 
smaller than for the single brane model. That is, naively, the answer from (i) and (ii) will 
be different than that from (iii). 

The resolution to this puzzle is already evident in the global structure analysis just 
discussed, where a zero of the metric corresponds to a particle horizon. There, we see that 
the metric used is geodesically incomplete, so we should not integrate beyond the particle 
horizon. In calculating Gn for the ith brane using the first approach, we should integrate 
over the 5th dimension only to the nearest particle horizon (i.e., the zero of the metric) 
on each of the brane. This also applies to the normalization of the trapped gravity mode 
(or any other mode) in the second approach. This way, the three approaches produce the 
same answer for Gn, as should be the case. Furthermore, Gn and the effective cosmological 



42 



constant is the same for all observers within the the same particle horizons, as pointed out 
earlier. 

The reasoning presented below may be summarized as follows. It takes infinite time for 
a signal emitted from a brane to reach the particle horizon, as seen by observers on the 
brane. The trapped gravity mode on the brane is treated as a perturbation of the metric in 
the "almost" conformally flat background metric |7j . The nearest particle horizons in the y 
coordinate map to z — ±00 in the "conformal metric" coordinate. Thus, the normalization 
of the gravity mode over the entire z coordinate only covers the region between particle 
horizons in the y coordinate, so the properties of branes separated by particle horizons are 
treated independently (wavefunctions and so Gn are independently normalized). 

For the negative tension brane in the Randall-Sundrum model, and for the visible brane 
in the compactified model discussed above, a naive determination using approach (iii) yields 
a negative G^, contradicting the above result. In Appendix |D|, we discuss the ambiguity 
and the problem arising in this approach. 

This above result also means that we cannot compare the mass scales between branes 
that are separated by particle horizons. In this case, the mass hierarchy problem in this 
brane world scenario is no worse than the standard hierarchy problem. We can still compare 
mass scales between branes that are not separated by particle horizons, and between mass 
scales measured by observers on the same (visible) brane. This implies that solutions to the 
mass hierarchy problem can be addressed only between branes that are not separated by 
particle horizons. Our discussion below follows closely the original discussion of Ref ||[7j . 



B. Green's Function and Newton's Constant 



Instead of bringing the metric into the conformally flat form ( |110| ), we can bring the 
metric (|36|) into an almost "conformally flat" form, where the y coordinate becomes the z 
coordinate given by: 

dz = ± dy = ± - ^ ^ - Vo ^ (126) 

f sinh [k (y - y )) H sinh [k (y - y )} ' 

Here we have in mind that H is very small, and all measurements involving Gn are at 
distances much smaller than 1/H. This above equation can be readily integrated, and the 
result is: 

k 

H (z + z ) = -lntanh- |y-y | ( 12 7) 

where zq is defined so that when y = we also have z = 0. We choose the "— " sign so 
that when y — > y we will have z — > 00. We observe that the space between the horizons 
( — Voi +l/o) is mapped into the entire real line in the new coordinate system. If the model 
contains many branes separated by horizons, each interval will be mapped into the entire 
real line, so we end up with a collection of disconnected spaces, each space containing one 
single brane. Using the above relationship between y and z we obtain the conformal factor 
of the metric: 

7^' '*<»-»>' = Wg (« + *)] • (128) 
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According to the result from [§,0] we obtain the following wavefunction for the trapped 
graviton mode: 



ij (z) = N 



H 



k sinh [H (z + Zq)] J 
where N is a constant that is determined by the normalization condition 

\^{z)\ 2 dz = 1. 



(129) 



(130) 



Calculating the integral we obtain 



2N 2 H 2 




1 + fc 2 



H 2 



■ 



1. 



;i3i) 



Expanding for small ff/A; we obtain: 



- I H 2 H 2 , /8P v 



(132) 



In the limit H — > we obtain iV = v/fc. 

In order to obtain the relationship between the 5-D and 4-D Newton constants, we first 
note that the relation between the physical perturbation and the trapped gravity wavefunc- 
tion is ||4| 



h 



^ 'physical 



(x,z)o^h^(x)A^/ A A(z)^ (z) 



(133) 



It turns out that ipo(z) oc A( d ~ 2 " 4 for this mode, so that the physical perturbation is simply 
h IJiU (x)A(z). To get the right scaling of G e ff on the brane, one also has to consider the 
source term carefully. For a point particle confined to the zth brane and at rest at x = xo, 
only the 00 component of the energy momentum tensor is nonzero: 



Tr 



00 



M 

A(zi 



■^(x-XoWz-zO, 



(134) 



where the brane is at z = Zi, A(z) = A(zi) and M is the particle mass measured on the 
brane relative to the physical metric, A(yi)hn U . The corresponding metric perturbation is 
(apart from possible exponentially small corrections p7t^5| ) 



2G b \jj { Zi )\ 2 M 
A(zi)\x. - x | 
2G^m 



Xo 



(135) 



where m = M^A(zi) is the mass of the particle in the effective 4-D theory (A similar 
analysis on this issue appeared while this paper was being completed p6| . ) 
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To see what the effective Newton constant is for brane-bound observers using the physical 
metric, h av = A(zi)h uu , let us first change coordinates to a system in which the unperturbed 
4-D metric is Minkowski on the brane. Then it is clear that the mass produces a perturbation 
which is identical to h Q0 , so all we have to do is re-express the solution found above in the new 
coordinate system. Since a coordinate separation r corresponds to a distance R = r{LA = 
rJA(zi), and m = M\J A(zi), we find that the perturbation is h 0Q = —2GNA(zi)M/R i.e. 
it is identical to Newton's law, but with a constant of gravitation G e // = G^A{zi). 

We would also like to know the corrections to Newton's force law that would be deduced 
on different branes. We find that, for two masses Mi and M 2 separated by a distance r(Lj) 
on the brane at Li, 

G N A(Li)MiM 2 , 2C 2 
r 2 {U) [ r 2 (L,) J 

where the masses Mi and M 2 and the distance r(Li) are measured with the induced metric 
7^. Here we also include the correction term due to gravity KK modes, where C 2 [ki, k i+ i) 
is a function of the bulk cosmological constants on the two sides of the Planck brane. We 
shall determine this function in the next section. When the the bulk cosmological constants 
on the two sides of the brane are equal, C 2 = k~ 2 . In terms of the rescaled metric 7^, 

rrii = <J A{Li)Mi and r(Li) = rJ A(Li), and relative to the induced metric there is a force 

This reproduces the familiar Newton's gravitational law as given by the 4-D effective the- 
ory fl59|). We see that the correction to Newton's force law on the visible brane becomes 
important at TeV scale distances. The Newton's force law in different metrics is discussed 
in terms of Feynman diagram in Appendix [E]. 

Consider the scenario in Figure 3. The above observation shows that the Newton's 
constant as seen by observers on each brane in Figure 3 yields the same relation ( |125| ), even 
though the ratio A(0)/A(L) can be exponentially large. The metric blow-ups beyond the 
particle horizons are unobservable. In fact, Gat as seen by observers on the two branes in 
the S l /Z 2 orbifold model will have the same value. By now, the reason is clear. In changing 
to the almost "conformally flat" metric, we have to determine separately for each brane, 
with each region between the horizons being mapped into real line, —00 < z < +00. Each 
region will contain a trapped graviton, and they do not "communicate" with each other 
across the particle horizon at y — yo since they belong to separate spaces. For an observer 
on the brane at y — 0, it takes infinite time for a light-like signal to travel from the brane 
to y : 

At = — {In [tanh k (y - y initia i)] - In [tanh k (y - y final)]} (138) 

where yinitiai = and At is clearly divergent when yfi na i —> yo- In terms of the 5-D 
parameters k 2 and qv, we see that A(0)/A(L) can be exponentially large. Since this huge 
factor is not measurable by observers on either brane, it has nothing to do with the hierarchy 
problem, that is, it neither solves nor worsens the standard hierarchy problem. 
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(136) 



IX. CORRECTIONS TO NEWTON'S GRAVITATIONAL LAW 



In the two brane model where both the hierarchy and the cosmological constant problems 
are solved, we note that the bulk cosmological constants on the two sides of the Planck 
brane are generically different. Here we consider this asymmetric situation, in particular 
the correction to Newton's gravitational law and the leading post-Newtonian effects. Let us 
summarize the key observation. The results are similar to that in the symmetric case, that 
is, the leading order post-Newtonian effects are the same as that from Einstein theory, and 
the correction to Newton's gravitational law is still r~ 2 as in the symmetric case 0, with 
a different coefficient. We shall calculate this coefficient below. Since we consider the case 
where the effective 4-D cosmological constant is exponentially small, it is safe to ignore it 
(i.e., set it to zero) in most applications. This simplifies the analysis. The analysis below 
follows closely that of Randall and Sundrum and that of Garriga and Tanaka f37H39 



Once we have the Green's function, we can also determine the correction to Newton's law 
for masses on the visible brane. Although the 5th dimension is compactified in the model, 
the gravity KK spectrum is still expected to be continuous because the model has a particle 
horizon. 

The Planck brane at y = is a surface layer separating two AdS spaces with different 
cosmological constants. We want to change the coordinates to Gaussian normal coordinates 
so that we can apply Israel's junction conditions. We expect the required change to be of 
the order of the perturbation induced by the mass located on the brane. Let us denote the 
new coordinates by: 

x a ± = x a ± + C±(x a ). (139) 

In general we will have to make different changes in coordinates for the two sides of the 
brane, so from now on we will suppress the ±. In the new coordinates the metric will be: 

ds 2 = g^dx^dx" + 2g fl5 dx fM dy + g 55 dy 2 . (140) 

The metric transforms as: 

9ab = 9 Pq g^g^; = 9ab + g aq i q b + g pb C% + O ((C) 2 ) ■ (141) 

We include now fluctuations around this background. We want to choose the functions 
£ a ( xb ) so that the perturbations h a b satisfy the RS gauge condition while the coordinates x a 
are Gaussian normal. Consequently we impose the conditions: 

h 5 5 = , = , ^5=0, g 55 = 1 

h 55 = 0, h ll5 = 0, = 0, h^ v = °- 

We transform now from the Gaussian normal coordinates to the Randall-Sundrum coordi- 
nates: 



ds 2 = 

-5 



g»u + gQi\\t*v) + dx^dx v + 2 \g^ + + g 55 ^ + dx^dy + 

^55 + 2^55^,5 + ^5 5 1 dy 2 = [g^ + h^] dx»dx u + 2 [g^ + h^} dx^dy + [g 55 + h 55 ] dy 2 . (142) 
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The condition h 55 = imposes £ 5 5 = so this function can depend only on the coordinates 
parallel to the brane, x p . The condition = imposes the condition: 



(143) 



since we keep only first order terms. This fixes the ^-dependence of the functions we 
have to impose continuity of these functions at y — 0, (0) = £^ (0) : 



2k 



(144) 



The functions £ A 
functions £ 5 (x p ) 



will be fixed by choosing the harmonic gauge on the brane, while the 
will be determined by the tracelessness of the RS perturbations. In the x a 



coordinates the background metric is given by: 



lvv (i - 2ke + ■ ■ 



(145) 



so we can obtain the relationship between the perturbations on the two coordinate systems: 



h 



2k 



V- 2 H-£ 5 -^ + W ( 146 ) 



Now the metric must be continuous at the brane. Expressing the metric in the Randall- 
Sundrum coordinates, the background metric is continuous at y = 0, so the fluctuations 
must also be continuous at y = 0. Using the fact that the functions £ M must be continuous 
at y — 0, and that £ 5 is independent of y, we obtain the following condition: 



k + e + (x p ) = k.e {x») . 



(147) 



We can understand this condition in terms of the "brane-bending" effect: the bending of the 
brane should be the same as an observer in the bulk approaches the brane from either side, 
y > or y < 0, so we can "glue" the two AdS spaces together along the brane. The relation 
( p.47 ) can also be understood in the following way. There are 5 physical polarizations for 
the 5-D graviton. For the massive KK modes, this is precisely what one needs for a spin-2 
particle. For the massless mode, the scalar component, if present, will contribute a Brans- 
Dicke-like interaction that is clearly ruled out. Fortunately, this scalar mode is gauged away 



by e 



Since there is only one scalar mode to be gauged away, Q. must be related. 



The jump condition at the brane, y = 0, can now be written as: 



Or, 



2k 2 



a 



(7^ + h 



:i48) 



where is the energy-momentum tensor and T its trace. Using the explicit form of the 
background metric and the fact that the bulk cosmological constants are determined to give 
an (almost) flat brane, the above equation can be simplified to the following form: 



\ \d r ,h 



y ltj [iu |0+ 



— d,-,h 



+ (k+ + k-) h^v = -k 



(149) 
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To obtain this equation we used the Csaki- Shir man solution for the brane tension 



401 



25 



2 2 

K a 



(150) 



The following steps are the same as in the paper |37 |: Combining the equation of motion 
for the metric perturbations and the jump condition, we obtain the Green function: 

{a- 2 n^ + d 2 - 4 (k\9 (y) + k 2 _6 (-y)) + 2{k + + 5 {y)}G R (x, x') = 5 5 (x - x') . (151) 

In the non-orbifolded case the Green function will have contributions from both "symmetric" 
and "antisymmetric" KK modes. In the orbifolded case the "antisymmetric" modes are 
projected out. The tracelessness of the Randall-Sundrum metric fluctuations will result in 
the following equation for the £± (x p ) functions: 



—k 2 T 
3 



□ < 4 ) 



k + e + + k_e 



2k 2 



k+k_T_ 
k+ + fc_ 3"' 



(152) 



We will use this result later when we calculate the corrections to Newton's law. The other 
four functions, £ M (x p ) provide the gauge freedom needed to restore the 4-D linearized Ein- 



stein equations. We will choose the usual harmonic gauge PT 

1 



8" h 







(153) 



so that the linearized Einstein equations for different components of the metric will decou- 
ple. Notice that the harmonic gauge condition is satisfied "up to KK corrections". The 
contribution from Q. (x p ) will change the numeric coefficient of T£ will change from | to 5 
so that the linearized Einstein equations become: 



□ 



(4) 



-2k 2 



T, 



2^ T 



(154) 



The set-up consists of a single positive tension brane placed at z = and the bulk is 
made of two slices of AdS^ with different cosmo logical constants, Aj. The background metric 
can be expressed as: 



ds 2 = e~ B ^ 



i] ab dx a dx b + dz i 



1 



r] ab dx a dx b + dz" 



(hz + 1, 

On each side of the brane, the RS graviton satisfies the usual Schroedinger-like equation 



(155) 



- d 2 z ip (z) + V (z) if> (z) = m 2 ij (z) 



(156) 



where: 



Viz) = —B'(zf--B"(z) 



(157) 



B(z) = 2\n(k s \z\ +k A z + l) 



(158) 
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where ks = + k + ) /2 and k A = (k + — /2, so that the equation for the wave function 
becomes: 



o2 / / \ , 15(k A + k s sgn(z)) 

-d z if) (z) + — —2 - Sk s S (z) 

4(k s \z\ + k A z + l) 



if) (z) = m 2 ip {z) ■ 



(159) 



In order to solve the above equation, we first make the change of variable, u — ki\z\ + 1 and 
then the change of function if) (u) = \/ux (u) so that the equation becomes: 



u d u x + ud u x + 



m 



u 



4 X = 0. 



(160) 



The m = mode has the following solution: 



ifj(z) 



z<0 

(fc_|z|+l)2 



N 

(k+\z\+l)V 



z>0 



161) 



The normalization constant N gives the 4-D Newton constant, 

2h+ k. 



G 



N 



Ge~ B ^ 4 \iP(0)\ 2 = GN 2 = G 



k, + k_ 



(162) 



In the limit k + = k- = k we obtain the RS relationship between G and G^: Gn = kG 

For the continuum mode, m^O, the solutions of Eq. (|159|) are linear combinations of 
the Bessel functions J 2 (jfru) an d Y 2 so the most general solution of Eq. (|159|) is: 



ifj(z) 



af_{z) + bg^(z) , z < 
c/+ (2) + d# + (z) , z > 



where: 



z| + 7— J 2 ( m (\z\ + ^- ) ) , g± (z 



k 



■^—Yo (m ( \z\ + -— 



fc. 



(163) 



(164) 



± , 



The coefficients a, 6, c and d are determined by the continuity condition at z = and the 
wavefunction normalizations. The calculation is straightforward but a little tedious. The 
result allows us to determine the correction to the Newton's gravitational law for the Planck 
brane: 



V(r) 



G^m\m2 



00 G N m e 



k 2k r 



-dm 



G^m\rri2 



1 q2 
1 + 2^2 



(165) 



where 



C 2 



1 



1 



4/4 1 - P + P 2 



10p(l -p + p 2 ) + 3(l + p 4 
3p(l + p 2 ) + 2p 2 



(166) 



where p = k + /k_. In the symmetric case (p = 1), C 2 = k 2 . For observers on the visible 
brane at L, 
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V(r) 



G jsimxTrii 



C 1 



2A(L)r- 



+ 



(167) 



where A(L) ~ 10 -30 is the warp factor. The physical fluctuations of the metric will be given 
by: = h^ + j^ k + Q_ + A;_£i where h nu is determined by V (r). Using the relationship 
Gn = 2k + k_G / {k + -\-kJ) we obtain the metric perturbations: 



h 



2G N M 



oo 



1 + 



2C 2 



3A(L)r 2 



hi 



2G N M 



1 



C 2 



3A(L)r 2 



8 



1.7 



(168) 



where L can be or Li These solutions differ from the ones obtained by expanding the 
Schwazschild solution to linear order by a coordinate transformation. 



X. DISCUSSIONS 

The viewpoint developed in this paper leads to a connection between the cosmological 
constant and mass hierarchy problems, resulting, rather generally, to relationships of the 
form (H/m pi anc k) 2 ~ {jnuiggs /mpianck) 2p ', with p > 1, when the expansion of the Universe 
is dominated by cosmological constant. Such a relationship is reminiscent of the so-called 
"large numbers" of cosmology, which have been variously either dismissed as numerology or 
regarded as physically significant clues to some fundamental connection between the micro- 
physical world and cosmology. Most explanations of how these numbers could arise from 
physics have tended to invoke variations on the theory of 4-D gravitation, such as advo- 
cated by Dirac |^2| or Brans and Dicke [ 43fl . However, Zeldovich [|4| proposed a different 



idea, that it is the cosmological constant, A e fj, not H that should appear in the large num- 
bers, and that they therefore tell us something fundamental about the relationship of A e ff 
to the rest of physics. The viewpoint developed here is intermediate, since it is based on 
multibrane solutions in 5-D gravitation, but only involves static spacetimes. In the context 
of such solutions, \n(m 2 Planck / m 2 Higgs ) depends on some combination of dimensionless sep- 
arations between branes, and \n(m Planck / A e //) depends on a different combination of the 
separations, and is always larger than \n(m 2 Planck j 'm 2 Higgs ) . The values of these dimensionless 
numbers depend on the particular arrangements of branes, which are not determined in these 
static models, and we do not yet understand what determines these arrangements funda- 
mentally. Moreover, even if we regard ^(rn 2 Planck /m 2 Higgs ) as observable, so one combination 
of brane displacements is knowable, the other combination, \n(m planck / A e ff), depends on a 
different combination, and is not fixed. It is clear that for dynamical spacetimes, we should 
expect some movement of the branes, resulting in variability of both ln(m Planck / A e f /) and 
^(rn 2 Planck /m 2 Higgs ). For some reason, which we do not address here, the branes must have 
been relatively stationary since before the epoch of cosmological nucleosynthesis, avoiding, 
in particular, changes in ^(m 2 Planck /m 2 Higgs ) that could alter the predicted light element 
abundances. 

In a 4-D spacetime theory, the cosmological constant includes contributions from micro- 
physics, which has many scales. For example, among other contributions, supersymmetry 
breaking around the electroweak scale will introduce a cosmological constant of the or- 
der of TeV 4 , which is many orders of magnitude bigger than the value observed, around 
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A e ff ~ (lO -12 ^) 4 . This is the cosmological constant problem. This problem will generi- 
cally persist if we consider a higher- dimensional theory and compactify it to 4-D, because 
the resulting effective 4-D theory will face the same fine-tuning problem. This problem may 
be avoided by considering higher dimensional theory with non-factorizable geometry, and 
the Randall-Sundrum model is the simplest example. In fact, the RS model has an (expo- 
nential) warp factor, which provides a solution to the well-known hierarchy problem. In Ref 
M, the RS model is extended to solve the cosmological constant problem. A key feature of 
the model is that the effective 4-D cosmological constant seen by an observer on the brane 
is not the vacuum energy density on the brane, but a non-linear combination of the brane 
tension and the bulk cosmological constants on the two sides of the brane. If other branes 
are relatively far from the visible brane, that non-linear combination becomes exponentially 
small, and the observer thinks he/she is seeing an exponentially small cosmological constant. 
For observers on the brane, it is natural to find an effective 4-D description. Then one may 
ask what happens to the argument for the need of fine-tuning in 4-D theory? The issue 
turns out to be quite interesting. Let us consider the single brane case first. If one insists 
on a 4-D description for observers on the brane, one ends up with a novel description via 
the AdS/CFT correspondence PlJlOll. Recall that, besides the trapped graviton, there is 



a continuous spectrum of Kaluza-Klein modes from 5-D gravity. Thus, in addition to the 
standard model fields (and whatever other fields that may exist) on the visible brane in an 
effective 4-D theory, an appropriate conformal field theory (CFT) must be included on the 
3-brane, to account for the effects of the gravity KK modes. That is, the CFT is equiva- 
lent/dual to the 5-D gravitational effects on the brane. This CFT interacts strongly and 
couples to other brane modes only via gravity. In fact, the classical correction to Newton's 
gravitational law in the 5-D picture becomes a quantum effect (i.e., the CFT contribution 
to the vacuum polarization of the graviton) in the 4-D picture [|45|J . Let us take the strongly 
interacting SU(N) N = 4 super Yang-Mills theory as the CFT. In the symmetric case, the 
bulk cosmological constant is related to the c- number of CFT 27t 4 k 2 A = 27 (N 2 — 1). 

Note that iV is discrete but not A. In this paper, we have explored the idea of adjusting 
the bulk cosmological constant in the context of static solutions of Einstein's equations, but 
it is reasonable to expect that similar adjustments occur in a dynamical spacetime. When 
branes move slowly, we would expect the bulk A to vary slowly and continuously, perhaps 
just making transitions along a sequence of static models like the ones constructed here. As 
A varies continuously, we have to consider another CFT when 27t 4 k 2 A = 27 (N 2 — 1) cannot 
be satisfied for integer N. In the uncompactified multibrane model, we expect the AdS/CFT 
correspondence to continue to be valid. Thus, as branes move, the theory explores the space 
of strongly interacting CFT theories, not just the parameter space of a given quantum field 
theory. That is, the effective 4-D theory changes as the brane separation changes. One 
should no longer apply here the argument for the need of fine-tuning inside a given 4-D 
theory to obtain a very small 4-D cosmological constant. 



XI. SUMMARY AND REMARKS 



In this paper we argue that, for parallel 3-branes that are relatively far apart, an expo- 
nentially small cosmological constant is quite generic. We present a two brane model which 
solves the hierarchy problem a la Randall-Sundrum and has an exponentially small cosmo- 
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logical constant for observers on the visible brane. Moreover, the cosmological constant, in 
Planck units, turns out to be much smaller than the Higgs mass in Planck units, naturally 
and generically. To be explicit, we write down the 4-D low energy effective action for the 
model. We consider this model as another motivation for the brane world. 

We see that the inequality (|J) is robust. This follows from the fact that both the visible 
brane and the Planck brane must be inside the same particle horizons. The cosmological 
constant is exponentially small roughly as a function of the distance of the Planck brane from 
the particle horizon while the hierarchy factor is exponentially small roughly as a function 
of the distance of the Planck brane from the visible brane. Since the visible brane must be 
between the Planck brane and the particle horizon, the inequality (|4j) follows. 

We discuss the physical implications of the particle horizons that may appear in the AdS 
bulk. In the process, we clarify some confusing issues in the literature. We also discuss a 
possible "phase transition" in the two brane model, which deserves more attention. We have 
concentrated exclusively on bulks that are purely AdS. Similar, simultaneous explanations 
for the cosmological constant and mass hierarchy problems may also be possible when the 
bulks are not pure AdS; this question also needs further study. 

There are interesting questions that remain to be addressed. Are the solutions stable 
under quantum effects ? We believe so. Both the hierarchy problem and the cosmological 
constant problems are solved by involving factors that are exponential functions of brane 
separations, so they become exponentially small as brane separation becomes large. This 
behavior resembles the Yukawa force, which is robust under quantum corrections. This 
leads us to believe that quantum corrections will not destroy the exponential behaviors in 
the brane world. 

It will be interesting to find a reason why the 3-branes like to stay parallel. The dynamics 
of these brane world models will be important. String/M theory realization is another 
important question. 
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APPENDIX A: THE STATIC SOLUTION TO THE BRANE EQUATION OF 

MOTION 



We can now check that the metric (fi2|) satisfies all the equations of motion derived from 
the action (|10|). The Einstein equations, the equation of motion of the field A ai ... ad _ 1 , and 
the equation for the induced metric have already been checked. Here we check the equation 
for the embedding coordinates, i.e., the brane equation of motion Eq. fllil) . In this appendix, 



we shall express the brane equation of motion in terms of the metric (|36|). There are three 
types of embedding coordinates, the spatial coordinates tangential to the brane, the time 
coordinate and the coordinate normal to the brane. We consider them separately. 
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The a = j = 1, 2, 3 case. 

This case is trivial, since all the terms in Eq. (14) are zero. The last term contains 
the field-strength tensor with two indices repeated: -P^ 1 ... M4 e Atl "' M4 and is zero due to 
the antisymmetry of the tensor. The first term can also be proven to be zero. The 
embedding coordinates are scalars with respect to the intrinsic coordinates of the 
brane, so we obtain: 



V M V„F 



7n| 



\l\ln\d£i 


-ir? 


\l\ln\ 




VlTnl 









(Al) 



since |7„| depends only on £° and Y n . We now prove that the second term of the 
equation is also zero: 



1 ab°a^ VpJl 7„ — L ab°a°l3 1n ~ 1 afi9 ~ cfl 
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since the metric is diagonal and does not depend on 







dx a 



(A2) 



The case a = 0. In this case the term containing the field strength tensor is also zero 
due to repeated indices and antisymmetry, but the first two terms are no longer zero. 
However, it will turn out that the two terms will cancel each other. The first term can 
be easily calculated: 
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(A3) 



In terms of the metric fl36|), this is equal to — ^ > . The second term can be calculated 
as in the previous case: 



F a/9 _ ± Qa 

1 apy — 2^ 
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(A4) 



which equals A( - L ; in the metric (|36|). The two terms cancel each other, so this equation 
is identically satisfied. 



The case a = 5. If a = 5 the term containing the field strength tensor is no longer 
zero, so we calculate this term first: 



F 5 r) Y ai ...P) V«4pMi-M4 



fJ.l---m.1l4 



e w-M4 = g^Fa^e^ = 4!F 50123 = 4!F 0123 5. 



(A5) 
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The equation of motion for the field A has already been solved, and the solution is: 



01— as _ e £ oi— 05 (A6) 



We can use this result to calculate the value of F01235: 

F01235 = F a ^g lai ■ ■ ■ g 5a5 = F 01235 g u ■ ■ ■ g 55 = 

9 F 01235 = -J\g~\ee 01235 = -eJ\g~\. (A7) 



The second term is: 

1 
2 



rld a X a dpX b ^ = Y%g^ = -g 5a 



dg a p dg aa dg a(3 



d£ a dtf d£ a 



1 pdgal = _lldg = _J_^J\g[ 

2 y dy 2gdy s fU dy ' 1 } 



Finally using the fact that JY) n \ = J\g (L n ) |, the equation becomes: 



dJ\g\ , — 

Vn—r, h e n ed\g\ = 0. (A9) 

9y v 

In the case where the derivative of g^ v has a jump across the brane, we should take the 
average of the values on the two sides. This yields Eq. (p9|) . Naively, we expect this 
equation to impose a non-trivial constraint between the brane tension and the brane 
charge, but in the static case we are interested in, this equation will be identically 
satisfied if the Einstein equation is satisfied. 

A similar result can be obtained for the brane equation of motion in the other models 
discussed in §|I[ 

APPENDIX B: DETAILS OF THE TWO BRANE MODEL 



Using the continuity and junction conditions ( 8IIJ84) , we can solve for each of the coth 
functions in (R3) in terms of ki,k2,qo,Qi, giving 

coth [A* (L - yi )) = " 2 \ % , coth [ kl (Li - yi )} = h_h±3L 

U2 i-2 „2 L.2 h.2 1 „2 

coth [k 2 {L X - y 2 )\ = fcl \ 91 , coth [k 2 (L 2 - y 2 )} = *> o \ + % . (Bl) 

2q 1 k 2 ' 2q k 2 

Following from Figure 2 and the subsequent discussion, the possible values of the various coth 
functions can be deduced. Using (|B~T|), these inequalities are constraints on the allowed values 
of ki and k 2 , given q and q%. Here coth [k% (L — yi)] < —1 and coth [k 2 (L 2 — y 2 )] > 1 always 
for positive tension branes, which imply k\ < (ki — q ) 2 and k\ < (k 2 — go) 2 respectively. 
The other inequalities depend on the behaviour of the scale factor in the bulks between the 
branes: 
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if there is a horizon between L\ and Lq only, 



coth [h (Lx - Vl )] > 1 
coth [k 2 (Lx - y 2 )] > 1 



k\ < (h - qx) 
k\ > (k 2 + q x ) 



(B2) 



if there is a horizon between L 2 and L\ only, 

coth [h (Li - yj)] < -1 = 
coth [k 2 (Li - y 2 )] < — 1 



fc| > (An + gi) 
&i < (h - qi) 



(B3) 



if there is a horizon in both bulks, 

coth [A* (Lx - y x )} > 1 =>• < 
coth [fc 2 (L x - y 2 )] < -1 => k\ < (k 2 



Qi) 

■Qx) 



(B4) 



The allowed regions of the "k-plane" can be further constrained by solving ( |BT|) for the 
brane separations in terms of kx and k 2 : 



L\ — L c 



2kx 



In 



Lx = — In 
2k 2 



(h + k 2 + qx) (kx - k 2 + gi) (fci + k 2 + q ) (h -k 2 + q ) 
(kx + k 2 - g ) (kx - k 2 - g ) (fci + A; 2 - gi) (fei - k 2 - qx) 
(k 2 + kx + qx) (k 2 - h + qx) (h + kx + g ) (k 2 - fci + g ) 
(fc 2 + /ci - gi) (k 2 -h- qx) (k 2 + kx - g ) (k 2 - kx - g )_ 



(B5) 



When there is no horizon between L and Lx, L 2 — Lx is positive is the region allowed by 
(|B3|), but the condition Lx — L > requires that k\ > kf + q qx. There is a reflection of this 
constraint about the line kx = k 2 in order for L 2 — Lx > when there is no horizon in the 
second bulk. Finally, we have imposed L 2 —Lx > Lx—L to distinguish the physically distinct 
configurations of the branes and horizons between them. The function (L 2 — Lx) — (Lx — L ), 



obtained from (|B5|) , can be evaluated numerically, revealing the regions in which this last 
condition is satisfied. Combining all constraints gives Figure 7. 

We find in Section that when there is a horizon in the second bulk, in the limit of 
large A/ x 2 , the hierarchy scale factor is fn 2 Higgs /m 2 Planck ~ exp[— 2Alx], and the cosmological 
constant scale factor is A e ff/rrip lanck ~ expj— (Alx + A/ 2 )] (with a horizon in the first bulk 
only, the roles of Alx an d AZ 2 will be reversed). We now prove that in our model, the 
hierarchy scale factor is always greater than that for the cosmological constant. The curve 
along which A/ 2 = Alx can be obtained from ( B5|) most simply by defining s = k 2 + kx and 
t = k 2 — kx- The condition AZ 2 = Alx becomes 

(t + qx)(t + q ) = (t-qx)(t-q ) 
(t-qx)(t-q ) (t + qx)(t + q ) 

which has t = or k 2 = kx as its solution for positive brane tensions. Above this line 
A/ 2 > Alx- Then for k 2 > kx, the only allowed region in Figure 7 is that for the single 
horizon in the second bulk; there A/ 2 > Alx, making the hierarchy scale larger. For k 2 < kx 
there are two allowed regions, one of which contains a horizon in both bulks, and hence 
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does not provide a solution to the hierarchy problem; the other region has a horizon in the 
first bulk only. In that region A/ 2 determines the hierarchy scale and All determines the 
cosmological constant scale, and since All > AZ 2 , the hierarchy scale factor is again the 
larger. 

A calculation similar to that which led to fl86|) reveals that, when there is a particle 
horizon in each of the two bulks, 

F 2 (Lx) ~ exp[-(Ah + Al 2 )\ (B6) 

with a complicated coefficient of proportionality. The Hubble constant in such a case will be 
greater than that with a single horizon because k 2 is less, giving a smaller A/ 2 ; this is evident 
from Figure 7. Hence the single horizon model is an energetically favorable configuration, if 
such a comparison is meaningful. 



APPENDIX C: DETAILS OF COORDINATE TRANSFORMATION 

In this appendix we describe the derivation of the coordinate transformation (|103|) - (|106| ). 
We start by scaling out all the dimensional constants by defining t = Ht, y = ki(y — yi), 
and x l = Hx l . The metric ( |101| ) can now be written as 

kids 2 = dtf - sinh 2 y di 2 + sinh 2 y e 2t d± 2 . (CI) 

If we focus attention on the first two terms in this line element, we see that the coordinate 
singularity at y = is analogous to the coordinate singularity of the two dimensional Rindler 
metric —y 2 dt 2 + dy 2 , and can be eliminated by the coordinate transformation 

y = tanh(y/2) cosht 

i = tanh(2//2) sinht (C2) 

The metric now takes the form kfds 2 = Q 2 ds 2 , where 

4 

n " = (i-f + py ^ 

and 

ds 2 = -di 2 + df + {i+y) 2 dx 2 . (C4) 

The coordinate singularity at y = has now been eliminated in the first two terms of the 
metric, but persists as the coordinate singularity i+y — of the last term in (|U4[). However, 
the metric (|C4| ) is just 5-D Minkowski spacetime 0. In double null coordinates u = t + y, 
v — t — y it can be written as 



5 We note that for the metric 

k 2 ds 2 = dy 2 - A(y)dP + ^(y)e 2 *dx 2 , 

if we define coordinates y = /(y)cosht and t = /(y) sinht and choose (f'/f) 2 = 1/A(y) (where 
prime denotes differentiation with respect to y), then kfds 2 = (f')~ 2 ds 2 , with ds 2 given by Eq. 
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ds 2 = —dudv + u 2 dic 2 . (C5) 

The transformation to Lorentzian coordinates for this flat metric can be obtained by solving 
for the exponential map from the tangent space at the point u — 1, v — 0, x l — to the full 
spacetime. This yields the coordinate transformation 

U = u 

V = V + MX 2 

X 1 = ux\ (C6) 
and the metric (|C4|) now takes the simple form —dUdV + dX 2 . If we now define coordinates 



T and Y by U = T + Y and V — T — Y, and combine together all the successive coordinate 
transformations of this appendix, we obtain the transformations ( |103|) - (|T06|) and the metric 
( |109| ) given in the body of the paper. 



APPENDIX D: G N AND A e// 

Here we explain why the determination of Gn from cosmology via the Hubble constant 
is sometimes misleading. Recall that H 2 = n%A e ff/3, where the 4-D gravitational coupling 
is given by Eq. fl56|), while the 4-D effective cosmological constant A e ff is given by Eq. (p8|). 

A eff = E A 2 (L t )a t - F(ki, L^, (Dl) 

so both k 2 n and A e ff are dependent on the positions Lj of the branes within particle horizons 
and the piecewise constant A, or equivalently, the fcj. Since a i: k? N) Lj and k t are related 
by Einstein's equation, a perturbation on one of the brane tensions, — ► oi + 5o \, requires 
a corresponding change in some of the other quantities. This back-reaction must be taken 
into account. Let us see when this back- reaction is expected to be small and when it may 
be sizable. Let us first consider the single positive tension brane model, where the brane is 
not charged. Here, H 2 = (q 2 — 4k 2 )/ '4, K 2 A e ///3 = q — 2k and k 2 n = 2kn 2 . We can vary 
the brane tension, or equivalently q by a small positive amount, q — > q + 5q. Keeping k 
constant, we find that dH 2 /dq = q/2 > and k 2 n does not change. This implies that the 
dA e ff/dq > 0. So we see that k 2 n = 35H 2 /SA e ff. As we shall see, a similar result will be 
obtained in a multibrane model if the metric is peaked at q brane. 

Now, let us consider a model where the hierarchy problem can be solved. Besides the 
visible brane, the model must involve at least another brane, where the metric factor A(y) 
is (exponentially) bigger. To be specific, let us consider a two brane model, where <j\ is the 
visible brane tension and <To is the Planck brane tension. 

Aeff = a + A 2 (Li)«7i - F(h, k 2 , L u L 2 ) (D2) 

where A{L\) is exponentially small. So the direct visible brane tension contribution to 
A e ff is exponentially small, and a small change in <j\ is totally negligible. However, as 
<7i — > <7i + Sax, either do, A(Li) or F(ki, k 2 , L%, L 2 ) will adjust accordingly, depending on 
the details of dynamics of the model. Since the influence of o\ on the other quantities are 
not necessarily exponentially suppressed, the induced change in A e jf can be large. It is easy 
to see scenarios where increasing a\ actually causes A e ff to decrease, or where A e ff still 
increases, but GV decreases in a way such that H 2 decreases. 
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APPENDIX E: NEWTON'S FORCE LAW 



Here we consider the Newton's force law between two point masses on the visible brane 
explicitly. An easy way is to compare it to the electric force law, since the electromagnetic 
field which is confined on the brane does not undergo a rescaling. 

Let us consider the electromagnetic field and two charged scalar fields confined on the 
visible brane at y = L: 

S (4) = / rf 4 x V^[i? (4) /24 - + £ ^Drf+D'h - A(L?M*%/2\. (El) 

3 

In solving the hierarchy problem, we expect A(L) ~ 10~ 30 . The electric charge e in = 



— ieA^ is of order unity. As before ||, we can rescale the scalar fields <pj = J A{L)<f>j and 



the masses rrij = ^A(L)Mj (J = 1,2) to obtain 

S (4) = / d*x^f[RW/2K% - -F, V F^ + £ \D^D^ - (E2) 

which is in the standard familiar form. Note that there is no rescaling for gauge field A^. As 
pointed out before, the effective action is independent of the way we define the warp factor, 
A(0) = 1, and A(L) = {m Higgs /m Planck f , or A(Q) = {m Planck /m Higgs f and A(L) = 1. 
The effective action is also invariant to any further changes in the coordinates. The 
difference in the coordinates parallel to the brane, ~ x b ) 2 ) wm represent the physical 

distance measured by an observer on the brane only if A(L) = 1, so we may use this 
convention. 

To warm up, let us first use the action ( |E2|) to calculate the Newton's force law by 
considering the one graviton exchange between the two scalar fields in the low energy ap- 
proximation with linearized gravity. In this large distance case, the energy-momentum tensor 
is T^ v ~ r]^ u J2j m |0| + .... The result is well-known: 

m\m\ 2 mim2 ,„ . 

— Crjv — 5 — , e — - — (Ldj 

where we also give the one-photon exchange case. Comparing to the usual electric force law 
e 2 /r 2 , we notice that we must factor out the m\m2 factor in the one-graviton exchange term 
to obtain the usual Newton's force law —G^mxm^jT^. This is the result in the (7^, <f>j, rrij) 
frame. Suppose we use the action ( [BTD instead, that is, in the (7^, <f>j, Mj) frame. The 
one-graviton exchange now gives 

_ ^ {A{LfMl){A{LfMD ^ 



r 



In the one-photon exchange case, we note that the coupling term ieA{L)A^cf) + d i j i (f) implies 
the electric charge in this case to be eA(L). Using the Fourier transform of A^Vjd^d^cf) — 
A(L) 2 M 2 (f) ~ 0, we see that 



; V /I(L)M 1 )( V /I(Z)M 



eA(L)Y ^ ' v - . (E5) 
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To obtain the electric force law e 2 /r 2 , we have to factor out a MiM 2 A(L) 3 factor, so the 
Newton's force law is given by 



A{L)M l M 2 mims 
" G n ~ 2 = -^n—}— (^b) 

so we obtain the Newton's force law of the effective theory, as expected. Since both forces 
are proportional to r~ 2 , a rescaling of r does not change the result. 

We can also compare the relative strengths of the electric and gravitational forces between 
two particles: 



Fq i t T \ M\M 2 



--G N A(L)-y-. (E7) 



Although the gravitational force is so much weaker than the electric force in the visible 
brane, the ratio is increased by a factor of about [pipi anc kj mm gg ^) 2 in the Planck brane. 
The variation of the relative strength of the forces as a function of the brane separation is 
analogous to the running of the couplings in ordinary field theory. 
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